PROBLEM 


Investigate  theoretically  the  dynamic  characteristics  of  underwater 
bodies  towed  ai  speeds  of  sfcoot  45  knots.  Specifically: 

1.  Devise  a  general  procedure  for  designing  towed  underwater  bodies  of 
optimum  dynamic  performance. 

2.  Present  a  practical  design  fulfilling  the  design  requirements. 

3.  Design  systems  for  automatic  control  of  both  longitudinal  and  lateral 
motion. 


RESULTS 

General  guidelines  for  the  design  of  an  underwater  towed  body  were  derived 
from  an  analysis  of  the  pertinent  hydrodynamics. 

2.  A  practical  tow  body  was  designed.  For  a  speed  of  45  knots,  the  maximum 
diameter  is  30  inches  and  the  main-body  length/diameter  ratio  is  2.7. 

3.  The  cable  tie-point  location  has  a  great  influence  on  the  body's  stability, 
and  it  should  be  forward  and  above  the  body’s  center  of  gravity. 

4.  The  hydrodynamic  coefficients  can  be  predetermined  by  means  of  a  mathe¬ 
matical  model  which  can  be  checked  on  analog  or  digital  computers. 

5.  Important  parameters  determining  the  dominant  characteristic  equation 
roots  for  longitudinal  and  lateral  motion  are  the  cable  tension  and  the  damping 
coefficients. 

6.  Tow  bodies  can  be  easily  designed  to  have  longitudinal  and  lateral  dynam¬ 
ic  stability. 

7.  Automatic  control  systems  are  desirable,  and  proposed  systems  are 
described  herein.  An  automatic  longitudinal  system  can  hold  the  tow  body  at  a 
commanded  depth  and  improve  the  stability  and  damping.  An  automatic  lateral 
system  can  improve  the  stability  and  damping,  coordinate  the  rudder  and  ailerons, 
and  provide  a  yaw-rate  directional  control. 

8.  The  attainable  depth  and  required  horsepower  depend  primarily  on  the 
towing  cable  and  on  the  body  drag.  Small-cross-section,  laminar-type  streamlined 
cables  are  required  to  achieve  low  drag  and  to  prevent  cavitation. 


RECOMMENDATIONS 

1.  Consider  the  use  of  the  body  designated  Body  A"  in  this  report,  as  a 
practical  tow  body, 

2.  Use  the  design  procedure  outlined  herein  to  design  additional  tow  bodies. 

3.  Consider  incorporation  of  the  automatic  control  systems  proposed  in  this 
report  in  tow  body  designs,  and  consider  development  of  the  necessary  hardware. 
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INTRODUCTION 

This  report  ewers  a  study  made  to  originate  (Da  general  design  procedure 
for  towed  underwater  bodies,  and  (2)  a  practical  body  and  cable  design.  The 
body  and  cable  should  be  such  as  to  enable  a  high-speed  surface  craft  to  drag 
a  specified  load  at  a  specified  depth.  If  a  sufficiently  large  horsepower  were 
available  at  the  surface  craft,  the  load  could  be  towed  at  the  specified  depth  and 
speed  by  any  stable  cable  of  adequate  length.  In  actual  practice,  it  becomes 
necessary  to  optimize  the  cable  and  body  design  so  that  (1)  the  power  require¬ 
ments  are  minimum,  (2)  the  ~tem  is  stable  so  that  the  body  stays  down  at  the 
prescribed  depth,  and  (3)  the  transient  characteristics  are  satisfactory. 

The  cable  design  has  to  be  closely  related  to  the  towed-body  design.  The 
cable  influences  the  body  dynamics,  changing  the  motion  modes  from  those  of  a 
self-propelled  body;  and  the  dynamics  are  also  significantly  dependent  on  the 
location  at  which  the  cable  is  attached  to  the  body.  Cable  requirements  are 
statistical  and  dynamic  stability,  minimum  drag,  no  cavitation  effects,  low 
vibration,  and  low  transmission  of  random  energy-  such  as  noise. 

The  dynamics  of  the  body  control  its  stability  and  transient  characteris¬ 
tics.  The  dynamics  can  be  treated  by  application  of  Newton's  second  law  of 
motion,  since  the  acceleration  of  the  body  is  resisted  by  hydrodynamic  forces, 
body  weight,  buoyancy,  and  cable  tension.  The  hydrodynamic  forces  can  be 
defined  in  terms  of  hydrodynamic  coefficients  that  depend  on  the  shape  and 
proportions  of  the  body.  These  coefficients  can  be  approximated  for  a  given 
body,  and  thus  be  used  to  form  the  basis  of  a  design  procedure.  The 
coefficients  and  the  effect  of  varying  them  are  examined  in  detail  in  this 
report. 

One  of  the  most  critical  design  requirements  is  minimum  drag.  The  ratio 
of  the  depression-force-needed-twhold-the-body-down  to  the  drag  has  to  be  made 
larger  as  the  depth  increases  and  as  the  towing  horsepower  decreases. 

The  maneuverability  of  the  surface  ship  with  its  towed  load  probably  is 
determined  principally  by  the  ship’s  characteristics.  In  this  report  we  examine 
the  desirability  of  automatic  controls,  and  propose  both  longitudinal  and  lateral 
control  systems.  The  longitudinal  system  tends  to  hold  the  towed  body  at 
commanded  depths.  The  lateral  system  coordinates  the  rudder  and  ailerons  and 
controls  the  yaw-rate  orientation.  These  systems  improve  the  stability  and 
damping  of  the  towed  body. 

Structural  vibrations  of  the  cable  and  body  are  mainly  randomly  excited. 

One  major  source  of  vibration  is  the  turbulent  boundary- layer  energy  around 
the  body  produced  by  its  high-speed  motion  through  the  water.  An  elastic 
cable,  and  a  compliant  surface  on  the  body,  can  act  as  filters  and  attenuate  this 
random  energy. 


DYNAMICS  ANALYSIS 


The  dynamical  properties  of  the  tow  body  can  be  explained  by  the  six 
equations  of  motion  that  express  Newton’s  second  law.  Three  of  these  equations 


equate  applied  force  to  the  rate  of  change  of  linear  momentum  with  time  in  the 
thee  space  directions.  The  other  three  equations  equate  applied  moment  to  the 
rate  of  change  of  angular  momentum  with  time  in  the  thee  space  directions. 
Although  these  equations  apply  to  a  body  moving  with  respect  to  fixed  inertial 
coordinates,  the  mathematical  solutions  can  be  simplified  by  a  transformation  to 
a  system  of  moving  coordinates  that  are  instantaneously  related  to  the  body.  In 
such  a  system,  the  moments  of  inertia  become  independent  of  time.  The  relation¬ 
ship  between  time  derivatives  referred  to  fixed  inertial  axes  and  those  referred  to 
moving  axes  is  given  by  tba  equation 


dV 

^fixed 


dV 

^moving 


+  coXV 


rotation 


where  V  is  s  linear  velocity 
u  is  an  angular  velocity 

The  moving  axes  consist  of  the  Cartesian  body  axes  instantaneously  fixed  in 
space  according  to  instantaneous  body  position.  The  derivative  term  is  position- 
fixed,  and  then  the  rotation  term  is  orientation-fixed.  The  time  derivative  of  the 
body  velocity  referred  to  the  moving  axes  appears  as  two  terms,  one  accounting 
for  translation  and  the  other  for  rotation.  The  equations  of  motion  when  written 
in  the  Sulerian  moving-coordinate  form  are  given  by  the  matrix  equations 
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where  F  is  the  applied  force  column  matrix 
M  is  the  applied  moment  column  matrix 
H  is  the  moment-of-momentum  matrix,  I co 
/  is  the  moment-of-inertia  matrix 
»  is  the  angular  velocity  of  the  body 


Six  Eulerian  coordinate  transformations  relate  the  space  orientation  of  the 
moving-body  set  of  axes  to  the  orientation  of  the  fixed  inertial  set  of  axes, 
through  three  angular  rotations.  For  our  purpose,  the  standard  aircraft  transfor¬ 
mations  of  yawing,  pitching,  and  rolling  have  been  chosen  (fig.  1).  These 
transformations  are  given  by  the  equation 
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FIRST.  POSITIVE  ROTATION  ABOUT  a  '  FOR  y  ANGLE 

NEXT.  POSITIVE  ROTATION  ABOUT  LINE  OF  NODES  FOR  0  ANGLE 

FINALLY.  POSITIVE  ROTATION  ABOUT  x  FOR  9  ANGLE 

Figure  1.  Coordinate  transformation  from  inertial  Cartesian  set  to  movin*-body  Cartesian  set. 


where  x,  y,  z  are  the  moving-body  set  of  coordinates 
x',y',z'  are  the  inertial  set  of  coordinates 

Y  is  the  Eulerian  angle  through  which  the  inertial  set  is  yawed 

0  is  the  Eulerian  angle  of  pitch  following  yawing 

9  is  the  Eulerian  angle  of  roll  following  pitching 

The  vector  angular  velocity  of  the  moving  body  is  the  sum  of  the  Eulerian 

angular  velocities.  The  yaw  rate  has  components  referred  to  the  moving  axes 
x,  y,  and  z  equal  to  -Y  sin  0.  Y  cos  0  sin  9  and  Y  cos  0  cos  9,  respectively. 
The  pitch  rate  has  components  referred  to  the  y  and  z  axes  and  equal  to  0  cos  9 
and  -0  sin  9,  respectively.  The  roll  rate  occurs  about  the  moving  x-axis  and 
has  only  an  x  component.  Thus,  the  vector  angular  velocity  w  has  x,  y,  and  z 
components  designated  p,  q,  and  r,  which  are,  respectively: 

p  =  $  _  4*  sin  0 

q  -  0  cos  +  ^  cos  0  sin  9 
•  • 
r  ~  Y  cos  0  cos  9  -  0  sin  9 

The  geometry  of  a  moving  submerged  body  can  be  described  by  the  use  of 
the  two  sets  of  Cartesian  axes.  The  inertial  set  is  chosen  so  that  the  x-  and 
z- axes  are  horizontal  and  the  positive  z-axis  points  downward.  The  moving-body 
set  is  chosen  so  that  the  positive  x-axis  points  longitudinally  toward  the  nose, 


the  positive  y-axis  to  the  right-hand  side  in  the  direction  of  motion,  and  the 
positive  z-axis  downward  as  in  a  right-handed  set  of  Cartesian  axes  (fig.  1). 

The  six  equations  of  motion  can  be  separated  into  two  sets  of  three 
equations  each,  one  set  defining  the  longitudinal  motion  and  the  other  the 
lateral  motion.  Longitudinal  motion  consists  of  motion  in  the  x-z  plane  and 
pitching  about  the  y-axis  without  coupling  energy  into  side  motion,  yawing 
motion,  or  roiling  motion.  Longitudinal  motion  is  defined  by  the  force  equations 
in  Fx  and  Fz,  and  the  moment  equation  about  the  y-axis  in  M.  The  yaw  and  roll 
rates,  and  the  side-slip  velocity,  are  assumed  to  be  zero.  The  lateral  motion 
is  defined  by  the  side-force  equation  in  F>t  and  the  yawing  and  roiling  moment 
equations  in  N  and  L.  The  pitching  moment  about  the  y-axis  is  assumed  to  be 
zero.  The  longitudinal  and  lateral  sets  of  equations  are: 


Longitudinal 
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is  the  sum  of  the  differential  forces  from  the  equilibrium  condition 
in  the  x-axis  direction 

is  the  sum  of  the  differential  forces  from  the  equilibrium  condition 
in  the  z-axis  direction 
is  the  body  mass 

is  the  body  moment  of  inertia  about  the  y-y  axis 

is  the  velocity  change  from  the  equilibrium  condition  in  the 

x-axis  direction 

is  the  velocity  change  from  the  equilibrium  condition  in  the 
z-axis  direction 

is  the  equilibrium  velocity  in  the  x-axis  direction 

is  the  angular  velocity  change  from  the  equilibrium  condition 

about  the  y-axis  of  the  body,  0 

is  the  total  differential  of  moments  about  the  y-y  axis 


Lateral 

2A  Fy  =  m(v  +  U0r) 
2AL  = Plxx~rjxz 
SA,V  -PJ„ 


where  2AFy  is  the  total  differential  of  side  forces  from  the  equilibrium 
condition 

2AL  is  the  total  differential  of  rolling  moments  from  the  equilibrium 
condition 

2aN  is  the  total  differential  of  yawing  moments  from  the  equilibrium 
condition 


p  is  the  roiling  anguiar-velocity  change  from  equilibrium 
r  is  the  yawing  angular-velocity  change  from  equilibrium 

Ixx  is  the  moment  of  inertia  about  the  x-x  body  axis 

/  is  the  moment  of  inertia  about  the  z-z  body  axis 

JX2  is  the  product  of  inertia  about  the  x  and  z  body  axes 

The  above  equations  have  been  written  for  small  perturbations  from  the  eqi 
rium  condition  of  straight-line,  constant-speed,  level  motion.  The  total 
differential  of  the  forces  and  moments  from  the  equilibrium  condition  balan 
the  change  in  momenta  from  their  steady-state  condition  (following  the  trer 
used  by  Bryan  and  Williams  in  their  1904  aircraft  stability  analysis).  The 
state  terms  which  balance  one  another  across  the  sides  of  the  equations  hi 
been  eliminated. 

Figure  2  show's  the  geometry  of  a  cable-towed  underwater  body  in  loi 
tudinal  motion  in  the  x-z  plane.  The  cable  tow  point  is  above  and  ahead  o 


body’s  center  of  gravity.  The  equilibrium  oqedkioa  i#  token  to  be  horiaoatol 
mot  km,  The  body’s  center  of  parity  and  center  of  bcoyancy  are  assumed  to  be 
at  the  sane  point  Thelongitudinal  equations  of  motion  written  in  terns  of  the 
hydrodynamic  force  and  moment  differentials  about  the  body’s  center  of  gravity 
become  for  the  *  component 

dF,  dFr  5F,,  BFX  5F  .  „ 

2AF  =  --*“  +  rrXw*  Tt*«’+-rf*  +  -x*-9  +  F=m(S-itq) 
du  dw  aw  3B  dv 


where  — *  u  is  the  force  differential  with  forward  velocitv 
du 

F  is  the  other  applied  force  due  to  the  cable  and  the  controls 
dF„ 

— •■w  is  the  force  differential  with  z-axis  velocity 
dw 

In  the  equation  above,  the  total  differential  is  the  sum  of  force  differen¬ 
tials  due  to  forward  and  downward  velocities,  force  differentials  due  to  attitude 
angle  and  angular  velocity,  and  force  differentials  due  to  the  cable  mid  tow-  body 
control  surface.  Each  term  is  in  units  of  force.  The  terms  are  divided  by  a 
quantity  with  the  dimensions  of  a  force  to  make  each  term  nondimensional.  All 
the  independent  variables  of  the  equation  themselves  are2divided  by  a  quantity 
of  the  same  dimensions  as  those  of  the  variable,  and  the  derivative  multiplier 
of  the  variable  is  multiplied  by  the  same  quantity.  Thus,  the  equation  can  be 
expressed  with  nondimensional  independent  variables  multiplied  by  nondiroen- 
sional  derivative  multipliers: 
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where  ~  is  the  change  in  angle  of  attack  la 


—  is  the  nondimensional  forward  velocity 
U 

S  is  the  lift  surface  area 


q  is  the  dynamic  pressure  defined  by  1/2  pU2,  p  being  the  fluid  mass 
per  unit  volume  and  l)  the  forward  velocity  (the  dimension  of  dynamic 
pressure  is  force  per  unit  area) 


c  is  a  representative  length,  the  lilt  surface  chord 


\J  is  the  disturbed  forward  velocity:  the  sum  of  the  equilibrium  velocity 
Uq.  and  the  velocity  variation  u 


The  liondimensioeid  iadepesoest  variable*  cm  be  written  Is  a  condensed  for*: 
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HYDRODYNAMIC  COEFFICIENTS 

The  noodimensionai  derivative  multipliers  are  known  as  stability  deriva¬ 
tives  or  hydrodynamic  coefficients  for  underwater  bodies.  They  are  written  with 
a  capital  C,  subscripts  indicating  the  particular  equation  of  motion  and  the 
particular  independent  variable.  These  coefficients  are  nearly  constant  over 
wide  ranges  of  speeds.  Their  values  depend  on  the  lengths  and  areas  used  in 
defining  them.  Sometimes,  the  maximum  cross-sectional  area  is  used,  but 
lengths  raised  to  powers  are  usually  employed  with  underwater  bodies.  The 
water  mass  density  does  not  always  feature  in  the  definition.  Here,  the  aircraft 
definition  with  lifting  surface  area  and  chord  is  adopted,  because  of  the  many- 
experimental  data  available  for  airplanes.  Coefficients  for  air  and  water  agree 
in  value,  when  defined  by  the  same  terms.  The  hydrodynamic  coefficients  iur 
the  x  equation  of  longitudinal  motion  are  as  follows: 
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Hydrodynamic  coefficients  can  be  ignored  if  they  represent  small  differen¬ 
tial  force  components  that  oppose  the  change  in  momentum.  The  '  ydrodynamic 
coefficient  for  the  x  component  forward-velocity-derivative  variations  is  small 
for  streamlined  bodies  with  a  large  enough  length-to-diameter  ratio.  This 
coefficient  is  called  the  virtual  mass  effect  and  represents  force  due  to  the 
change  in  momentum  of  the  displaced  water  mass. 

The  hydrodynamic  coefficients  can  be  estimated  with  reasonable 
accuracy,  and  they  can  be  used  as  the  basis  of  the  design  of  an  underwater 
lowed  body.  Here,  the  nondimensional  coefficients  are  given  per  radian  angle, 
also  a  nondimensional  quantity.  The  longitudinal  coefficients  whose  magnitudes 
are  large  enough  to  affect  the  dynamics  significantly  are  listed  in  table  1,  along 
with  their  definitions,  simplified  formulas  for  evaluating  them,  and  typical 
numerical  values  for  aircraft: 


TABLE  1.  LONGITUDINAL  HYDRODYNAMIC  COEFFICIENTS 
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TABLE  2.  (Continued) 


Coefficient 

Definition 

; 

I 

Formula 

Aircraft 
Typical  Value 

C-r 

i  2 u  as 

-V  „  s,JiSacL„ 

-0.1 

Sqb  b  dr 

4  ■  '“n  S  \b)  da 

1 

Sq  C0 

qt  $tv 

dat  q  S 

-0.5 

c *, 

1  dF 

y  a 

TSf  dS r 

c 

0.2 

CnSf 

1  dNa 

Sqb  dSr 

c* 

"P  dsr 

-0.1 

1  3La 

Sqb  dSr 

CL  Cnbr 

4  Cm 

0.01 

cma 

l  Mg 

Sqb  <96  0 

function  of  AR,  aileron  location, 
and  chord  ratio  to  wing  chord 

0.6 

c *. 

Sqb  dba 

- c, .  £», 

8 

-0.015 

Symbols  are  as  follows: 


b  is  the  lift  surface  span 
r  is  the  dihedral  angle  in  degrees 

y  is  the  distance  from  body  center  along  the  wings  to  the  center  of  pressure 
of  each  wing 

Sp  is  the  area  of  a  wing  with  a  dihedral  angle 
is  the  area  of  the  vertical  tail  surface 
t(v  is  the  moment  arm  from  the  vertical  tail  quarter-chord  to  the  body’s 
center  of  gravity 

CN  is  the  side  force  coefficient 

q  is  an  efficiency  factor 

F  is  the  side  force  due  to  the  control  surfaces 
Na  is  the  yawing  moment  due  to  the  control  surfaces 

La  is  the  rolling  moment  due  to  the  control  surfaces 

6r  is  the  rudder  deflection;  positive  rudder  causes  positive  yawing 
6q  is  the  aileron  deflection;  a  positive  aileron  deflection  causes  a  positive 
roiling  moment 


Approximate  equations  for  the  Interal  hydrodynamic  coefficients  cnn  b« 
formed  from  the  approximate  formulas  given  in  table  2  if  the  values  of  the  nearly 
constant  quantities  ar  taken  as: 
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di  q.  dCN  dC, 

—  =*  0.5,  r,  =  1,  =  1,—^  =  2,  — ^  =2 

aa  q  dat  a  a 

On  the  basis  of  these  values,  the  lateral  hydrodynamic  coefficients  become 
(approximately): 


-r 

~  57.3 

C(p  =  -0.5  ( AR  =  5.  elliptic) 

C«,= 

r 

-  2  $tv  ^ tv 

C .u 

C  -  L 

c  - 

S  b 

np  jg 

°nr 

Cy?> 

_  O  $tv 

Cy6r  ~~°-5Cyp 

cn6jr= 

ct6r 

..  CL  C»*r 

Cg  6  =  0  to  about  0.6 

cnfi  = 

4  C„r 

n6a 

£l 

4 


■u.o 

-Cl  c«L 


a  8  Ci 


£p 


LONGITUDINAL  EQUATIONS  OF  MOTION 

The  longitudinal  equations  are  defined  by  the  force  equations  in  Fx  and  Fz, 
and  the  moment  equation  about  the  y-axis  in  M.  The  force  equations  are: 


IA  Fx  =  m(u  +  wq) 

Z&F;  =  m(w-U0q) 

For  the  towed  underwater  body,  force  differentials  with  the  body  attitude 
angle  arise  because  of  body  weight  and  buoyancy  and  because  of  the  tow  cable. 
Figure  3  illustrates  the  resolution  of  the  weight  and  the  buoyancy  forces  in  the 
body  axes.  The  center  of  gravity  (c.g.)  and  the  center  of  buoyancy  are  assumed 
to  be  coincident.  The  x  and  z  force  components  are  then  given  by 


Fgx  ~  ~<ms  ~  B^sin  0 


Fgz  =  (mg  -  B) cos  0 


where  mg  is  the  body  weight 
B  is  the  body  buoyancy 

The  changes  in  these  forces  with  attitude  angle  are  given  by  the  derivatives: 


T5 


Figure  3-  Weight  nnd  buoyancy  resolution. 


3F. 
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—  =  -(mg  -  B >  cos  0 


dZ 

36 


—  =  -(mg  -  B)  sin  0 


The  tow  cable  produces  differential  forces  in  the  x  and  z  axes: 

AF xT  =  T  cos  0e  +  T0  sin  0C0 

AF^T  =~Tsin  0C  +  T0cos  0p0 

where  T0  is  the  undisturbed  cable  tension 
T  is  the  differential  cable  tension 
0C  is  the  undisturbed  cable  angle 
0  is  the  change  in  attitude  angle 

A  differential  moment  is  produced  about  the  y-y  axis  due  to  the  changes  in 
attitude  angle  and  cable  tension: 


AM  =  -T  cos  6cdz  +  T  sin  6cdx  -  TQ  sin  6cdzd  -  T0  cos  6cdx  0 

where  dx  is  the  distance  the  tie  point  is  ahead  of  the  body’s  c.g. 
dz  is  the  distance  the  tie  point  is  above  the  c.g. 

The  Laplace  transformations  of  the  small  perturbation  equations  of  motion 
can  be  written  in  terms  of  hydrodynamic  coefficients  that  are  the  force  and 
moment  derivatives  of  changes  in  angle  of  attack,  changes  in  attitude  angle, 


and  changes  in  forward  velocity.  Their  first-order  time  derivatives  become  for  a 
cable-towed  underwater  body  with  elevator  control: 


/mUs  \ 

c")'“ 


x  force  equation 

l-c  n  \  W-B 
+fo77c,&s~C*a)a  +-r-fcos8ie  +j 


(l  G  s_Io  sinA 

\2  U  xq  3q  c) 


T  1  dF*Se  e 

= —  COS0.+ - £  6. 

Sq  c  Sq  «?8„  e 


•-*;  ,_r  •  v® 

"  r  .. 

*  '  *,  #  >  vH 


2  force  equation 


I mU  c  (  \  n-mU  c  1 

-c»'“+  (|sj  'iuc-js  -  c»y  “  +U  C^J s  ~  sr  r°C089' 

i  \  -r  .  i  . 

+  —  W-B  sin  0  I  6  =  — -  sin  6  +  - - —  6 

Sq  I  Sq  c  Sq  d  5  c 


y  moment  equation 


-C  'ul  —  C  *s  -  C_„ 
mu  \  2{J  roa 


y  \Sqc  2  u 


Cms  +  — —  sin  0_cfe  + 
m<?  Sqc  c 


s^c°8vfe)1 


1  /  \  1  <5M6 

— j-T  cos  0„(fe  +  Tsin0  </x)  + - 2  8,, 

qc\  c  /  S?cd6, 


These  complete  longitudinal  equations  of  motion  can  be  simplified  if  the 
variations  with  forward-velocity  fluctuations  are  not  significant,  as  may  be  the 
case  in  constant-velocity  towing.  The  motion  can  then  be  described  by  the 
z  force  and  the  y  moment  equations,  with  'a  and  0  as  the  independent  variables: 


z  force  equation 


mV  c  \  f—mU  c  \  1  1 

- C  .  s-C  '<*  +  - -  —  C,„  s-—  T,.cos0„  +  —  W-B  sin0  0 

Sq  2V  za)  10  \  Sq  2 V  zq)  Sq  0  r  Sq  c 

-T  1  dFza 

=  —  sin0„  + - £  8  , 

Sq  c  Sq  dS  e 


17 


y  moment  equation 


LONGITUDINAL  EQUATIONS  FOR  SPECIFIC  TOWED  SODIES 


For  this  report,  three  specific  underwater  bodies  were  quantitatively 
studied.  A  design  parameter  basic  to  all  three  bodies  is  a  30-inch  maximum- 
diameter  cross  section  Below  we  detail  the  analysis  of  a  body  designated 
Body  A  (fig.  4), using  both  the  complete  longitudinal  equations  and  the  short 
equations  (for  insignificant  forward-velocity  fluctuations).  For  the  analysis, 
the  body  is  assumed  to  be  towed  at  45  knots  in  a  straight-line,  level  path.  The 
numerical  parameters  for  Body  A  are: 

U  =76  ft/sec  S  =  6  ft2  m  =60  slugs 

lyy  =  400  slug-fl2  Tq  =  14,000  lb  l(  =  6  ft 

x'  =6.8  inches  dz  -  1.167  ft  dx  =0.167  ft 

St  =2  ft2  9C  =  81.85  degrees  W-B  =  500  lb 

q  =  5776  lb/ft2 

The  hydrodynamic  coefficients  for  longitudinal  motion  can  be  obtained 
from  the  approximate  formulas  given  in  this  report.  Based  on  the  body  param¬ 
eters,  the  coefficients  become: 


CXU  =  -CD  =-0.12 

X' 

C  =4-  ■ 

ma  c 

e,  s, 

-1.5  -1  =  -0.7333 

c  S 

C,„  =  -CL  =  0.8 

c.j  =  -s(- 

?A2  S. 

J.  _L  =  -36.0 
c)  S 

C,a  =  CL  =  -0.4 

-4 

'.\2  S. 

1  _L  =-72.0 
c)  S 

dC. 

C  =-— -^  =-4-0 

*a  do. 

C„6  =->'5 

e 

C/i 

II 

L 

© 

i.  s. 

c„  =  -e7-|  =  -'20 

C'6c  =F7 

Cm6p  =-0.5 
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B.  SIDE  VIEW 


5  WING  SECTIONS  MODIFIED  NACA  FOUR  DIGIT 

I  0020-65 

j  5  INDICATES  MAXIMUM  AT  0.5  CHORD 

J  6  INDICATES  LEADING  EDGE  RADIUS 

00  INDICATES  SYMMETRY 
20  INDICATES  THICKNESS  0.2  CHORD 


MAIN  BODY  DTMB 
(LANDWEBER-GERTLER)  FORM 


C„  -  0.55 


m 


-=0.5 


J_ 

d 


p  n  R*t 

,o.V„0.45  r,.— -0. 


d  2 


Figure  4.  Form  of  Bodies  A.  A',  and  A”  (dimensions  in  inches). 


and  x‘  is  toe  distance  by  which  the  lift-surface  quarteMhord  is  ahead  of 
the  c.g. 

qt  is  the  dynamic  pressure  at  the  tail 

8  is  the  angle  of  deflection  of  the  elevator  (taken  positive  for  downward 
deflections) 

SM is  the  static  margin  (or  the  chord  ratio)  by  which  the  c.g.  is  ahead  of 
the  neutral  point 

Approximate  design  equations  can  be  derived  by  using  values  near  the 
following  for  some  of  the  quantities  that  do  not  vary  greatly  from  design  to  design: 


<k 

—  =  0.5,  K  -  1.0, 

da 


0.5, 


=  3.0, 


dCf 

da 


=  4.0 


Based  on  the  preceding  approximations,  the  longitudinal  hydrodynamic  coeffi¬ 
cients  Lecome: 


zq 


t.  S. 
.6-1  -i 
c  S 


C.i  = 


e.  s. 

-S-  — 
c  S 


Cxu  =  -2  CD 


mq 


='6(t) 


2s, 

s 


cmi  =  -31 


2  c 
t\ 


cl  S 


cxa=cL 


x'  i.  S, 

C»*=47-L5it 


c,  =  -- 


d  a 


CzB =  ~o~  Cm 6 
e  e 


l,  S, 

Cm6  =  -1.5—  rr  C,„  ~  -2  C, 


c  S 


Thus,  most  of  the  longitudinal  hydrodynamic  coeiTicients  are  functions  of 
the  moment  arm  from  the  quarter-chord  of  the  horizontal  tail  area  to  the  body’s 
center  of,  gravity  and  of  the  ratio  of  the  size  of  the  horizontal  tail  area  to  the 
lifting-surface  area.  The  magnitudes  of  the  approximate  coefficients  agree  fairly 
well  with  measured  values.  The  approximate  design  formulas  are  based  on 
assumed  values  such  as  0.5  for  the  slope  of  the  downwash  angle  with  angle 
of  attack,  dt  /da.  This  value  can  be  diminished  by  raising  the  tail. 

Below  are  given  the  complete  longitudinal  equations  of  motion  for  Body  A, 
with  numerical  values  of  the  coefficients  inserted: 


(0.13158s  -  0.12)  'u  +  (0.4)  ’a  +  (-0.29405)0 

=  4.0830  x  10-6  T 

(-0.8)  1  u  +  (0.17105s  i  4.000)  'a  +  (-0.05263s -0.057162)  0 

=  -2.8564  xlO"5T  -0.5  6 p 
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(0)  'u+  (0.23684s  +  0.73333) '«  +  (0.011542s2  +  0.47368s  +  0.38943)9 


=  -36 


e 


The  solution  of  these  equations  of  motion  for  the  homogeneous  case  of  no 
right-hand  applied  terms  gives  transients  of  the  same  exponential  form  for  all  of 
the  independent  variables  'a,  u,  and  0: 

S,/  „  So*  „  So*  „  s.  i 

u  =  Ae  1  +  Be  2  +  Ce  3  +  I)e  4 
where  A,  B,  C,  D  are  constants. 

The  quantities  sj,  s%,  S3,  »-4  are  the  roots  of  the  characteristic  equation, 
which  is  the  polynomial  equation  in  s  obtained  from  the  determinant  of  the 
coefficients  of  the  independent  variables  of  the  homogeneous  simultaneous 
equations.  The  value  of  the  determinant  is  zero  if  the  transient  solution  is  to 
exist.  The  roots  of  the  characteristic  equation  must  be  negative  or  have  negative 
real  parts  for  the  transient  solution  to  be  stable.  Stability  can  be  indicated  by 
the  Routh  Hurwitz  method  without  solving  for  the  roots.  Solution  for  the  roots 
directly  indicates  the  stability  of  the  solution.  Graeffe’s  method  is  an  accurate 
and  straightforward  means  of  solving  polynomial  equations.  It  is  based  on  raising 
the  polynomial  to  higher  powers  so  that  the  roots  separate. 

The  characteristic  equation  for  the  complete  longitudinal  equations  of 
motion  for  Body  A  is  a  quartic: 

s4  +  71.649  s°  +  1096.8  s  2  +  2496.2  s  +  1821.1  =  0 

The  roots  of  this  equation  obtained  by  Graeffe’s  method  are  a  pair  of  negative 
real  roots  and  a  pair  of  complex  roots  with  negative  real  parts: 

Sj  =-51.14 
s2  =-17.96 

s.j  =  -1.28  -  j'0.63 
s4  =  -1.28  +  j'0.63 

The  complex  roots  represent  a  dani|«d  sinusoid.  The  complex  pair  can  be 
written  as  a  quadratic  equation  in  terms  of  the  natural  frequency,  «,,,  and  the 
damping  ratio,  £ : 


s2  +  2  s  +  w„2  =  0 

For  Body  A,  the  natural  frequency,  coy,,  is  1.408  radians  per  second  mid  the 
damping  ratio,  £ ,  is  0.9055.  The  complex  terms  damp  to  half-amplitude  in  the 
time 

0.693 

T  =  -  seconds 
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The  real  root  transient  terms  represent  pure  exponentials  and  damp  to  half- 
amplitude  in  the  following  time: 


<P  — 


0.693 
root  value 


seconds 


The  time  to  damp  to  half-magnitude  for  the  complete  longitudinal  equations  for 
Body  A  is  determined  by  the  smallest  of  the  two  real  root  times  <0.0136  and 
0.0388  second)  and  the  complex  root  time  <0.544  i  ->nd). 

The  short  solution  of  the  longitudinal  equations  ignores  the  a  velocity 
fluctuations  and  is  determined  by  the  z  force  equation  and  the  y  moment  equation. 
The  characteristic  equation  is  a  cubic  and  becomes,  for  Body  A. 


s3  +  70.738  s2  +  1018. 12  s  +  7G9.7M  =  0 


The  roots  of  the  cubic,  obtained  by  Graeffe‘s  method,  are 
s,  =-51.18 
s2  =  -18-82 
s3  =  -0.80 

The  transient  solution  is  again  stable  but  consists  of  three  pure  exponential 
terms.  The  time  to  damp  to  half-amplitude  for  the  three  roots  becomes  0.0135, 
0.0368.  and  0.865  second.  The  values  of  the  first  two  are  nearly  those  of  the  two 
large  loots  of  the  complete  solution.  These  two  are  related  to  the  short  period 
times  of  aircraft.  The  third  root,  resulting  from  the  cable,  is  the  dominant,  one 
and  determines  the  body  damping.  A  study  of  the  characteristic  equation  of  the 
short  solution  shows  that  the  dominant  or  smallest  root  is  approximately  given  by 


T  .  T 

—  sin  9  dz  +  — -cos9  dx 
Sqc  c  Sac  c 


c 

2V  Cmq 


The  dominant  cable  root  can  be  increased  (damping  time  reduced)  by  increasing 
the  cable  tension  (depression  force).  The  cable  tie  point  should  be  in  the  same 
direction  as  the  depression  force-drag  ratio  of  the  body.  The  iow  point  should 
lie  ahead  and  above  the  center  of  gravity  of  ttie  body.  The  pitch-damping  hydro- 
dynamic  coefficient  should  not  be  too  large.  Increasing  the  tail  length  increases 
the  damping  time  of  the  dominant  root. 

For  Body  A,  the  time  to  damp  to  half-amplitude  is  larger  with  the  short 
solution  than  with  the  complete  solution.  Introducing  the  x  equation  of  motion 
makes  a  complex  root  of  the  dominant  short-solution  root.  An  approximate  value 
of  the  time  to  damp  to  half-amplitude  can  be  obtained  from 


c 


9.693  x  2  x  «y  CitfnQc 
2U 


Sq 


-~C (J-^sine  dz  +  ~costf.dxVc 
Sq  za\Sqr  e  Sqc  '  )  2 


C„  —  c_„  +  cc 

St  2C  2QQ  ZO  22/ 


2U 


mq 


This  formula  agrees  with  the  short-solution  approximation  in  that  it  shows  the 
desirability  cf  decreasing  the  pitch-damping  hydrodynamic  coefficient  and 
increasing  the  depression  force  in  order  to  decrease  the  time  to  damp  to  half¬ 
amplitude.  The  damping  time  is  also  decreased  by  increasing  the  drag  and  lift 
coefficients  of  the  body. 


LAPLACE  TRANSFORMS  OF  LONGITUDINAL  EQUATIONS 


The  Laplace  transform  equations  of  motion  under  applied  forces  are 
simultaneous  linear  algebraic  equations  that  can  be  expressed  in  terms  of  the 
independent,  variables  by  matrix  inversion.  They  can  be  algebraically  solved  by 
Cramer’s  rule: 


l C]T 

[cl 


f 


where  {/}  is  the  column  matrix  of  the  independent  variables  («,  ’a,  0 ) 

[C]r  is  the  transpose  of  the  matrix  of  cofactors  of  the  hydrodynamic 
coefficients 

|c|  is  the  determinant  of  the  hydrodynamic  coefficients 
{f}  is  the  column  matrix  of  the  applied  forces 

Transfer  functions  are  the  ratios  between  the  independent  variables  ar.d  the 
applied-force  variables  obtained  from  the  simultaneous  transform  equations. 
These  s-plane  variations  of  the  independent  variables  'u,  'a,  and  0  with 
changes  in  the  elevator  angle  and  the  cable  tension  are  given  for  Body  A  by 
the  following  equations: 
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In  the  transfer  functions  above,  the  tension  change  T  is  in  pounds,  the 
velocity  fluctuation  u  is  nondimensional,  and  the  angles  are  in  radians.  The 
Body  A  transfer  functions  for  the  complete  longitudinal  equations  are  shown  in 
Appendix  A,  and  those  for  Body  A'  (with  the  opposite  dihedral  angle)  are  shown 
in  Appendix  B. 

The  short-solution  transfer  functions  for  Body  A  are  also  shown  in 
Appendix  A.  The  s-plane  variations  according  to  the  short  solutions  are: 
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Examination  of  the  complete-solution  and  the  short-solution  transfer 
functions  indicates  that  there  are  two  basic  types  for  each  solution: 


I 

I 


Complete  Solution 

S2  +  OjS  +  Oq 

ts  +  y)  (s  +5.1  (fs  +  al2  +p2] 

s3+  a0 s  +  a,s  +  a0 
(s  +  y)  (s  +6 )  [ts  +  a)2  +  |32j 

Short  Solution 


s  +  a0 
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fs+  y)  (s+8Hs+M 


The  Laplace  transform  of  the  independent  variables  is  equal  to  the  product 
of  the  transfer  function  and  the  transform  of  the  applied  force.  Since  the  trans¬ 
forms  of  the  impulse  functions  equal  constants,  the  transforms  of  the  independent 
variables  are  of  the  same  form  as  the  transfer  funct  s  for  impulse-type  applied 
forces.  For  unit-step-function  applied  forces,  the  transforms  of  the  independent 
variables  become  the  product  of  the  transfer  functions  and  the  transform  of  the 
unit  step.  1  s. 

The  time  solution  of  the  independent  variables  can  be  obtained  by  using 
Laplace  transform  pairs.  For  the  short  solution,  the  angle-of-attack  and  the 
attitude-angle  variations  with  either  elevator  angle  or  cable  tension  impulse  force 
are  of  the  form 

Ac’*1  +  Be’6'  +  Ce'M 


where  A,  II,  C  are  constants  depending  on  the  transfer  function  constants 
y,  6,  h  are  constants  (roots  of  the  characteristic  equation) 

The  solution  consists  of  three  transient  terms  that  decay  to  zero  at  rates 
determined  by  the  roots  of  the  characteristic  equation.  The  time-solution  form 
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depends  primarily  on  the  denominator  of  the  transforms,  the  numerator  (which  is 
at  least  one  order  smaller  than  the  denominator)  only  influencing  the  constants 
A,  B,  and  C.  The  time  solution  for  a  step-function  applied  force  is  still  of  the 
same  form  but  adds  a  constant  term.  The  steady-state  response  is  zero  for 
impulses  and  constant  for  step  functions  as  can  also  be  seen  from  the  equation 
for  the  steady-state  response: 

fit)  =  sF  (s)f(s) 

t  —  oo  S  0 

where  F(s)  is  the  transfer  function 

f(s)  is  the  transform  of  the  applied  force 

The  above  equation  is  zero  for  impulse  functions  and  constant  for  step-function 
applied  forces. 

For  the  complete  longitudinal  equations,  the  time  solutions  for  the 
independent  variables  (u,  'a  ,  and  6)  for  impulse  applied  forces  are  of  the  form 

Ae~Yt  +  Be'6*  +Ce'al  cos  (p/+  q>) 

where  A,  B,  C  are  constants  depending  on  transfer  function  constants 
<p  is  a  constant 

a  is  the  damping  term  in  the  complex  root,  £wfl 
p  is  the  natural  frequency  term  in  the  complex  roots 
y,  6  are  roots  of  the  characteristic  equation 

The  solution  for  impulse  forces  consists  of  damped  exponentials  and  a  damped 
sinusoid.  The  solution  for  step-function  forces  is  of  the  same  form  but  adds  a 
constant  term.  The  steady-state  response  is  zero  for  impulses  and  constant  for 
step-function  inputs.  Thus,  the  towed  body’s  longitudinal  response  to  impulses 
consists  of  transients  that  decay  to  zero  while  the  body  follows  step-function 
amplitudes. 


LATERAL  EQUATIONS  OF  MOTION 

The  lateral  small  perturbation  equations  of  motion  can  be  written  in  terms 
of  the  three  independent  angle  variables  of  side  slip,  yaw,  and  roll: 


2AFy.  =  m('0(p  +  T> 

<y  force  equation) 

SAL  =?/„-?/„ 

<x  moment  equation) 

2AV  =  ?/„-*$/„ 

(y  moment  equation) 

with  r  =  p (  0 

and  ( ’  =  L() 
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where  v  is  the  side-slip  velocity 
Uq  is  the  undisturbed  velocity 
i!  is  the  disturbed  velocity  along  the  x-x  axis 
P  is  the  angle  of  side  slip 

/  is  the  moment  of  inertia  about  the  x-x  axis 

/  is  the  moment  of  inertia  about  the  z-z  axis 

/  is  the  product  of  inertia  about  the  x  and  z  axes 

Control  forces  and  moments  caused  by  both  ailerons  and  a  rudder  are  assumed. 
The  aileron  control  is  assumed  to  be  a  differential  system,  the  motion  of  the 
aileron  on  the  right  wing  opposing  that  on  the  left  wing.  The  aileron  deflection, 
6a,  is  one-half  the  sum  of  the  motion  of  the  up  and  down  ailerons.  Positive 
rudder  and  aileron  deflections  cause  positive  yawing  and  rolling  moments.  The 
significant  hydrodynamic  forces  and  moments  are  functions  of  the  independent 
angle  variables  and  their  first  derivatives.  The  body  weight  and  buoyancy,  and 
the  cable  tension,  produce  differential  forces  and  moments  due  to  the  changes 
in  angles  of  roll  and  yaw.  These  become  hydrodynamic  coefficients  when  they 
are  made  nondimensional: 


A  F..  -(mg-Blai 

>  W-ll 

Af\vT  =  <~T0  Sin0c>9 
AF,.t  Sf-T0  cos0r)T 
A Lt  =  (~T0  sinQ  cdz)<p 

A /VT  =  f-T0  cos0  cdxlT 

Figures  5  mid  6  show  the  geometry. 

With  hydrodynamic  coefficients  C},„  Cyr.  and  Cya  assumed  zero,  the 
Laplace  transforms  of  the  equations  of  motion  written  in  terms  of  the  hydrodynamic 
coefficients  and  the  control  forces  and  moments  become: 


y  Force 
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BODY  DIRECTION  CHANGES  DUE  TO  CHANGING  THE 
DIRECTION  OF  MOTION  (TURNING)  AND  DUE  TO  SIDE 
SLIPPING.  SO  THAT  THE  BODY  AXIS  IS  NOT  IN  THE 
DIRECTION  OF  MOTION 


INITIAL  MOTION 

TURN  ANGLE 


SIDE  SLIPPING  IN  DIRECTION 
OPPOSING  TURNING 


FINAL  MOTION 


p  ANGLE  MEASURES  THE 
CHANGE  IN  BODY  DIRECTION 
DUE  TO  SLIP;  POSITIVE 
CLOCKWISE  MEASURED  FROM 
FINAL  BODY  DIRECTION  TO 
DIRECTION  OF  MOTION 


y  ANGLE  MEASURES  CHANGE 
OF  BODY  DIRECTION  WITH 
TURNING  AND  SIDE  SLIPPING; 


ANGLE  OF  YAW  IS  NEGATIVE 
OF  ANGLE  OF  SIDE  SLIP 


Future  5.  Motion  in  the  x-y  lnternl  plane 


GRAVITY  FORCE  DUE  TO  ROLL 


CABLE  MOMENT  DUE  TO  ROLL 


CABLE  TENSION  DUE  TO  YAW 


CABLE  MOMENT  DUE  TO  YAW 


Figure  (!.  I  .literal  forces  anti  moments  due  to  roll  and  yaw 
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LATERAL  EQUATIONS  FOR  SPECIFIC  TOWED  BODIES 

The  numerical  analysis  of  the  lateral  equations  of  motion  is  detailed  below 
for  Body  A,  assuming  that  the  body  is  towed  at  45  knots  in  a  straight,  level  path. 
The  numerical  parameters  of  Body  A  are: 


U  =76  ft/sec 

S  =6  ft2 

m 

=  60  slugs 

l22  =  400  slug'll2 

T0  =  14.000  lb 

l(D 

=  6  ft 

1  xx  -  45  slug-ft2 

dz  =  1.167  ft 

dx 

=  0.167  ft 

lX2  -  -100  siug-ft2 

c  =  81 .85  degrees  W'-B 

=  500  lb 

Stu  =  1.5  ft2 

q  =  5776  lb/ ft2 

b 

=  6  ft 

Based  on  these  parameters,  the  hydrodynamic  coefficients  are: 

Op  .;£j  -  -0.0524 

Cnr  =  -4 

1 

tl 

04 _ 

4.0 

C(.p  =n AR,  Shape)  =-0.5 

c,v  = 

£l 

4 

-0.1 

c  P 

r  -  o  tl'  <»'  -nt 

"P  S  b 

ctdr~ 

C"6r  _ 

4  r„r 

-0.025 

c, 

C„„  =  -  -ii  =  0.025 

»P  16 

c„y 

0  5  ^’/ip  =  “ 

-0.25 

f,r4  =-0-5 

c< 

-<5.  CV 

8  9, 

-0.06 

C  5  =  0.5  C  A  =0.25 

CfV 

- 

0.6 

The  lateral  equations  of  motion  for  Body  A  with  the  numerical  values  of  the 
coefficients  inserted  are  ns  follows: 

y  Force 


0.294059  +  (0.13158s  +  0.057102)4'  +  <0.1315hs  +  0.5>p  =  0.258, 


1  *  ' 


r<; .  ■  i 


x  Moment 


(0.00021641s2  +  0.019737s  +  0.077781)9  +  (0.00048092s2  +  0.0039474s)? 


+  0.0524P  =  -0.0256r+  0.66  Q 


z  Moment 


(0.00048092s2  +  0.00098684s>9  +  (0.0019237s2  +  0.039474s  +  0.0015910)?-  0.5p 

=  -0.256r  -0.065 a 

The  solution  of  the  homogeneous  equations  with  no  applied  right-hand 
forces  gives  transient  exponentials  and  a  damped  sinusoidal  term  for  each  of  the 
independent  variables  <p,  ?,  and  p.  The  roots  of  the  characteristic  equation 
obtained  from  the  nonzero  determinant  of  the  multipliers  of  the  independent 
variables  determine  the  exponential  constants  of  the  transients.  The  towed-body 
characteristic  equation  is  a  quintic.  Typical  airplane  characteristic  equations 
differ  from  those  of  the  towed  body  in  that  they  contain  no  constant  term.  The 
zero  root  indicates  that  the  plane  is  insensitive  to  the  choice  of  heading,  there 
being  no  moment  tending  to  return  it  to  the  original  heading.  The  towed  body’s 
path  is  determined  by  the  cable.  The  characteristic  equation  for  Body  A  is 


s5  +  247.48  s4  +  6689.3  s3  +  85,111  s2  +  246,333  s  +  92,743  =  0 


Solution  of  this  quintic  characteristic  equation  by  Graeffe’s  method 
indicates  stability,  since  there  are  three  negative  real  roots  mid  a  pair  of  com¬ 
plex  conjugate  roots  with  negative  real  parts.  The  roots  for  Body  A  are: 


s  j  =  -0.442 
s2  =  -3.34 

s3  =  -13.04  +  j  10.867 
s4  =  -13.04  -  j  10.867 
s5  =  -217.63 


The  transient  solution  for  lateral  motion  is  given  for  each  independent  variable 
by  an  equation  of  the  form: 

9  =  Ac'0-442'  +  Be3-341  +  Ce’217-63'  +De‘13-04t  cos  16.97/ 


The  times  to  damp  to  half-amplitude  for  each  of  the  terms  are  1.57,  0.2075, 
0.0531,  mid  0.00318  seconds.  The  two  small  roots  dominate  the  damping  of  the 
transients.  Examination  of  the  equations  of  motion  shows  that  the  largest 
damping  time  (smallest  root)  is  approximated  by 
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t  =  0.693 


+  C 


nU 


Cos8crnp 


62  mV 

mtv  nr  +  ~&T 


cosSj. 


The  time  to  damp  to  half-amplitude  is  decreased  for  the  smallest  root  by 
increasing  the  cable  tension  and  decreasing  the  yaw  damping  coefficient, 
Cnr.  The  time  to  damp  to  half-amplitude  for  the  second  smallest  root  is 
approximated  by 


— —  sine  JzC„r 
Sqb  c  nr  Sq 


mV  b  b  b  mV 

2V  +  CfpW  *~Sq  °lp  2V 


b 

f-c 


"P 


t=0.693  • 


Sqb 


b  mV  Tft 

Tnsin6  dz —  C  +C  -  - - —  sin6  ,(h 

0  c  2V  nr  Sq  Sqb  c 


This  next-to-smallest  root  resembles  the  spiral  divergence  root  for  aircraft.  The 
cable  tension  replaces  the  rolling  moment  coefficients.  The  damping  time  is 
decreased  by  increasing  the  depression  force  and  decreasing  the  roll  damping 
coefficient,  Cfp.  The  dihedral  angle  somewhat  affects  the  size  of  the  smallest 
root,  but  generally  does  not  affect  the  two  smallest  roots  as  much  as  in  the  case 
of  aircraft. 

If  straight-line,  level  motion  is  assumed  with  the  side-slip  angle  equal  to 
the  yaw  angle,  and  the  roll  angle  and  its  derivatives  being  zero,  the  yawing 
moment  equation  becomes 


-'n5. 


,2_.  * 


Sqb 


2V 


Cs  +  (C,.o  +  — —  cos  dx) 

nr  ,,p  Sqb  c 


For  Body  A.  the  denominator  is 


0.0019237  s2  +  0.03947-1  s  +  0.5015910  =  0 


The  roots  of  this  quadratic  are 


s,  =  -10.26  +;1 2.47 
s2  =  -10.26  -  /12.47 

The  damping  coefficient  (,  is  0.635,  and  the  i  atural  frequency  oo>,  is  16.15.  For 
the  complete  equations  of  motion,  the  complex-pair  roots  give  £  =  0,768  and 
=  16.97.  The  complex  roots  represent  the  “Dutch  Roll.” 

If  the  side-slip  and  the  yaw  angles  are  assumed  to  be  zero  with  no  applied 
rudder  forces,  the  rolling  moment  equation  can  also  lie  directly  solved: 
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lxx  9  ^  j.  r0  ■  J 

—  s*  -  ■ —  s  +  — sin  dz 
Sqb  2  U  Sqb  c 


For  Body  A,  the  denominator  of  the  transfer  function  is 


0.00021641  s2  +  0.019737  s  +  0.077781 

The  roots  of  this  equation  are  -4,13  and  -87.08.  The  transfer  function  equation 
indicates  that,  for  the  towed  body,  the  steady-state  roll  angle  is  zero  with  an 
impulse-type  aileron  deflection.  The  steady-state  roll  angle  becomes  a  constant 
for  a  step-function  aileron  deflection.  The  ailerons  act  like  elevators  in  that,  for 
a  steady-state  response,  the  deflections  must  be  held. 


LAPLACE  TRANSFORMS  OF  LATERAL  EQUATIONS 

The  Laplace  transform  lateral  equations  of  motion  can  be  solved  in  terms 
of  any  of  the  independent  variables  <p,  T.or  p.  with  either  rudder  or  aileron 
applied  forces,  by  Cramer’s  rule.  The  transfer  functions  which  are  transform 
ratios  between  the  output  mid  the  input  variables  are  numerically  evaluated  for 
Body  A  below: 
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The  solution  of  the  transforms  in  time  can  l>e  rejiresented  by  exponentials  and  a 
damped  sinusoid  for  impulse  inputs.  The  solution  adds  a  constant  to  these 
transients  for  a  step-function  input.  The  steady-state  response  to  an  impulse- 
type  input  is  zero  for  all  the  above  transfer  functions.  The  steady-state 
response  to  a  step-function  input  is  a  constant  for  all  the  transfer  functions. 

The  lateral  motion  analysis  above  was  based  on  Body  A  with  a  dihedral 
coefficient  of  -0.0524.  Consider  another  body.  Body  A’,  with  a  dihedral 
coefficient  of  +  0.0524  but  with  the  same  values  as  for  Body  A  for  the  other 
coefficients.  Based  on  this  sign  reversal,  the  roots  of  the  qusntic  characteristic 
equation  of  lateral  motion  become 

s  |  =  -0.353 
s2  =  -4.24 
«3  =  -12.42+;  11.69 
s4  =  -12.42-;'  11.69 
s5  =  -218.05 

For  these  roots,  the  times  to  damp  to  half-amplitude  are  1.96.  0.1634,  0.05580. 
and  0.003178  seconds,  compared  to  1.57  .  0.2075,  0.05314.  and  0.003184  seconds 
for  Body  A  with  the  negative  dihedral  coefficient  The  dominant  cable  root  is 
slightly  smaller  in  magnitude  while  the  other  two  real  roots  axe  slightly  larger. 
The  complex  Dutch  Roll  mode  has  a  slightly  smaller  real  part  but  a  slightly 
larger  absolute  magnitude  than  for  Body  A  Based  on  the  damping  time  of  the 
dominant  cable  mode,  the  body  with  n  negative  value  of  (Yp  appears  more 
desirable.  (In  the  control-system  analysis  that  follows,  it  is  indicated  that  the 
body  with  the  positive  dihedral  coefficient  may  be  preferable.) 


**wM*rv t 


The  transfer  functions  giving  the  variations  of  the  independent  variables 
and  £  in  Laplace  transforms.  in  terms  of  the  applied  rudder  and  aileron 
inputs,  arc  for  Body  A*: 
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CONTROL  SYSTEMS  GENERAL 

An  automatic  control  system  is  desirable  if  it  enhances  the  accuracy  of  an 
action  an  j  reduces  the  amount  of  needed  human  attention.  The  preceding 
analysis  of  the  towed  body  indicates  stability  for  both  longitudinal  and  lateral 
motion.  An  impulse  elevator  deflection  in  longitudinal  motion,  and  impulse 
rudder  or  aileron  deflections  in  lateral  motion,  all  produce  responses  that  return 
the  body  to  its  initial  reference  position  with  time  in  a  stable  manner.  The 
towed  body  follows  step-function  inputs  in  both  longitudinal  and  lateral  motion 
with  an  error  in  the  steady  state.  Automatic  control  systems  should  be  installed 
if  it  is  desired  to  tow  the  body  in  particular  motion  modes  and  if  it  is  necessary 
to  improve  the  dynamic  responses.  Some  advantages  that  can  be  obtained  with 
closed-loop  feedback  control  systems  are  listed  below: 

1.  Less  attention  is  required. 

2.  The  transient  response  can  be  improved,  so  that  the  body  returns  to  the 
reference  condition  after  a  disturbance  without  excessive  oscillations.  Transient 
characteristics,  such  as  the  time  to  rise  to  maximum  response  after  a  step  input, 
the  overshoot  above  the  steady-state  value,  and  the  time  required  to  settle  to  the 
steady-state  value,  can  be  adjusted.  Optimizing  [trocedures,  such  as  minimizing 
the  integral  of  the  product  of  time  and  the  absolute  value  of  error,  can  be  used. 

3.  The  response  to  specific  inputs  can  be  made  to  occur  without  error  or 
with  only  a  smali  error,  in  the  steady  state. 

4.  The  design  can  be  such  as  to  afford  good  stability  characteristics. 

The  stability  should  not  he  affected  by  small  changes  in  amplifier  gains.  After 
a  disturbance,  the  system  should  return  to  the  reference  condition  quickly. 

5.  The  control  ratio  between  the  output  and  the  input  can  lx;  made  less 
susceptible  to  undesired  disturbances.  The  control  system  can  be  designed  to 
adapt  itself  to  changing  conditions. 

6.  The  control  system  can  be  designed  such  that  the  response  to  random 
statistical-type  disturbances  is  minimum.  A  wide  frequency  bandwidth  in  the 
system  increases  the  random  disturbance  power,  but  allows  more  accurate 
following  of  rapidly  varyi  ig  input  signals  and  can  improve  the  system’s  stability. 
The  bandwidth  should  be  large  enough  to  puss  the  desired  input  frequencies,  but 
should  attenuate  the  higher  frequencies.  The  system’s  resonant  frequencies 
should  not  be  near  external  resonant  frequencies. 

In  a  single-loop,  feedback  control  system,  part  or  all  of  the  output  is  fed 
back  to  the  input.  (It  can  be  fed  back  in  any  phase  relationship  to  the  input; 
negative  feedback  occurs  when  the  feedback  signal  subtracts  from  the  input; 
positive  feedback  when  the  feedback  signal  adds  to  the  input.)  The  system  error 
Els)  is  the  difference  between  the  reference  input  and  the  signal  fed  back,  in 
negative  feedback;  and  is  the  sum  of  these  two  in  positive  feedback.  The  control 
ratio  between  the  Luplace  transforms  of  the  output  and  the  input  is  related  in 
the  transfer  functions: 


Cl  si  _  Gls) 
RIs)  ”  1  +  Gls)H(s) 


where  Cls)  is  the  transfer  function  of  the  output 
RIs)  is  the  transfer  function  of  the  input 


G(s)  is  the  forward  transfer  function  for  the  circuit  connecting  the  output 
to  the  system  error,  E(s ) 

H(s )  is  the  feedback-circuit  transfer  function 

For  negative  feedback,  the  algebraic  sign  of  the  G(s)H(s )  term  is  positive  for 
real  positive  G(s)  and  H(s)  values.  In  this  case,  the  control  ratio  when  H(s) 
equals  unity  varies  from  0  to  1  as  the  forward  transfer  function  G(s)  varies  from 
0  toco.  For  positive  feedback  with  G(s )  and  H(s)  positive  and  real  numbers,  the 
control  ratio  with  ll(s)  equal  to  unity  varies  from  0  to  «  as  G(s)  varies  from 
0  to  1,  mid  from  -ooto  — 1  as  G(s)  varies  from  1  to  ».  The  product  G(s)H(s) 
is  known  as  the  open-loop  transfer  function.  When  this  is  large  compared  to 
unity,  the  magnitude  of  the  ratio  between  the  output  and  the  input  becomes  equal 
to  1  /ll(s). 

The  system  error,  E(s),  for  negative  feedback  becomes  equal  to 

E(s)  _  1 

R(s)  ~  1  +  G(s)H(s) 


The  algebraic  signs  of  the  transfer  functions  G(s)  and  H(s)  are  assumed  to  be 
positive  in  the  above  equation. 

The  stability  of  closed-loop  control  systems  is  determined  by  the  denomina¬ 
tor,  1  +  G(s)H(s),  in  these  equations.  Instability  occurs  when  the  open-loop 
transfer  function,  G(s)IKs),  is  of  unity  magnitude  and  has  an  180-degree  phase. 
The  open-loop  transfer  function  can  usually  be  written  for  linear  systems  as  a 
ratio  between  factorable  polynomials: 


G(s)H(s) 


K(s  +  alls  +  b)(s  +  c)(s2  +  2£  t  +  cq,^  .  .  . 

sn(s  +  k)(s  +  70s  +  ml/s2  +  2£2co„2  +  ’  ’  * 


where  a,b,c  are  open-loop  zeros 


k,l,m  are  open-loop  poles 


K  is  the  gain  constant 

wf|  is  the  natural  frequency  part  of  a  complex  pair 
£  is  the  damping  ratio  part  of  a  complex  pair 


The  form  of  the  open-loop  transfer  function  indicates  the  complexity  of  the  control 
system.  The  denominator  must  be  at  least  one  degree  higher  in  s  than  the  numer¬ 
ator,  for  practical  systems.  When  the  denominator  does  not  exceed  the  numerator 
by  more  than  two  degrees  in  s,  the  phase  does  not  generally  equal  180  degrees  for 
finite  frequencies  unless  there  is  an  s 2  factor  in  the  denominator.  A  factor,  s", 
shifts  the  phase  n  times  90  degrees.  Factors,  s  +  a,  approach  90-degree  shifts  as 
s  =  jo}  increases  in  value.  Quadratic  factors  approach  180-degree  phase  shifts  as 
w  increases.  Open-loop  transfer  functions  with  the  denominator  exceeding  the 
numerator  by  one  or  two  degrees  in  s  usually  indicate  inherently  stable  control 
systems  unless  there  are  two  open-loop  poles  at  s  =  0.  When  the  denominator 
exceeds  the  numerator  by  more  than  two  degrees  in  s,  the  control  system  still 
is  usually  stable  over  wide  ranges  of  gain  values,  K.  As  the  gain  setting  is 
increased  from  zero  (the  open-loop  value),  the  closed-loop  [roles  change  (move 
along  the  root  loci)  and,  eventually,  instability  occurs.  It  is  desirable  that  the 
closed-loop  [roles  (roots  of  the  characteristic  equation,  1  +  Gtslflfsi  =  Q)  be 


37 


sufficiently  far  to  the  left  in  the  complex  plane  to  insure  rapid  transient  decay 
(large  negative  roots  and  large  negative  real  parts  of  the  complex  roots  of  the 
characteristic  equation).  Stated  in  another  way,  inherent  stability  occurs  if  the 
magnitude  of  the  open-loop  transfer  function  decreases  at  the  most  20  dB  per 
decade  in  frequency  or  40  dB  per  decade  with  the  sn  factor  in  the  denominator 
having  n  equal  to  either  0  or  1.  Each  degree  in  s  in  the  numerator  causes  the 
open-loop  transfer  function  to  increase  by  20  dB  per  decade,  while  each  degree 
in  s  in  the  denominator  causes  it  to  decrease  by  20  dB  per  decade.  The  phase 
of  the  open-loop  transfer  function  is  always  less  than  180  degrees  for  inherent 
stability. 

The  steady-state  error  of  the  closed-loop  systems  depends  on  the  form  of 
the  input  functions  in  addition  to  the  system’s  transfer  functions.  Control 
systems  are  classified  into  types  by  the  power  of  the  sn  factor  of  the  denominator 
of  the  open-loop  transfer  functions.  Type  0  has  n  =  0,  Type  1  has  n  =  1,  and  so 
on.  The  final  value  Laplace  theorem  can  be  used  to  determine  the  steady-state 
error.  The  theorem  states  that  if  the  Laplace  transform  of  f(t)  is  sF(s),  mid  if 
sF(s)  is  analytic  on  the  imaginary  axis  and  the  right-half  plane,  then  the  final 
value  of  the  time  function  is  given  by 

lim  f(l)  =  Jim  sF(s) 
t  ~+  »•  s-*-  0 


On  the  basis  of  the  final  value  theorem,  the  steady-state-system  following 
error  becomes 


e(t)  -  eft)  =  lim 
ss  t -*  oo  s  — *-0 


sR(s) 

T  +  G(s)H(s) 


For  a  Type  0  control  and  a  step  function  input,  the  steady-state  following  error 
is  a  constant.  For  a  Type  l  control  and  a  step  function  input,  the  output  follows 
the  input  without  error  in  the  steady  state. 

Control  systems  are  usually  modified  from  a  preliminary  design  to  improve 
the  stability,  the  transient  response,  and  the  following  error.  Compensation 
consists  of  changing  the  forward  transfer  function,  the  feedback  transfer  function, 
or  both  by  adding  phase  lead  or  phase  lag  networks.  Networks  added  may 
include  90-degree  phase  leads  such  as  differentiation  circuits.  90-degree  phase 
lags  such  as  integration  circuits,  and  RC  filters  which  are  integration  and 
differentiation  circuits  in  themselves. 

The  effects  on  the  control  systems  of  vurious  compensation  methods  are 
examined  below  in  a  general  manner.  High  damping  usually  indicates  better 
stability  characteristics.  Transient  damping  is  indicated  by  the  time  required 
to  damp  to  half-amplitude  for  either  simple  roots  or  complex  root  pairs.  The 
small  roots  (dominant  roots)  determine  the  least  damped  components.  More 
stability  is  indicated  by  larger  negative  real  loots  or  larger  negative  real  parts 
of  complex  root  pairs.  The  usual  effects  of  compensation  methods  are: 

1.  Integration  in  feedback  path.  Adds  u  closed-loop  pole  and  a  0.1  pole  at 
zero.  The  steady-state  error  is  less  since  the  system  is  one  type  higher. 
The  damping  and  the  stability  effects  are  usually  not  large  since  the  effects 
on  the  dominunt  roots  are  generally  small.  The  effects  can  be  large  for 
large  dominunt  roots  of  the  characteristic  equation. 

2.  Integration  in  forward  transfer  [unction  Adds  a  closed-loop  pole.  The 
system  is  one  type  higher  mid  follows  n  more  complicated  input  function 
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without  error.  Damping  and  stability  effects  are  usually  smalt  for 
small  roots. 

3.  Differentiation  in  the  feedback  path.  Removes  a  closed-loop  pole  and  adds 
an  open-loop  zero  at  zero.  The  system  is  one  type  lower  and  has  a  poorer 
steady-state  condition.  The  damping  and  stability  effects  are  usually 
small  for  small  roots. 

4.  Differentiation  in  forward  transfer  function.  Removes  a  closed-loop  pole. 
The  system  is  one  type  lower,  and  less  complicated  input  functions  can  be 
followed  without  error.  The  damping  and  stability  effects  are  usually 
small  for  small  roots. 

5.  Tachometer  feedback  H  +  b  s  feedback  loop).  The  steady-state-error 
constant  is  higher,  and  thus  the  steady-state  error  is  higher.  The  damping 
and  stability  are  both  improved. 

6.  High-pass  filler  added  in  series  with  the  differentiation  in  tachometer 
feedback.  The  steady-state  error  is  improved  over  tachometer  feedback. 

The  damping  and  stability  are  both  improved. 

7.  Lead  compensation  in  forward  transfer  function. 

Tt>T 2-  Adds  a  closed-loop  pole  and  a  closed-loop  zero.  The  steady- 
state  error  is  the  same  as  without  the  lead  compensation.  The  damping 
and  the  stability  are  both  improved. 

8.  Lag  compensation  in  the  forward  transfer  function 

T,<  T2-  Adds  a  closed-loop  pole  and  a  closed-loop  zero.  The  steady- 
state  error  is  not  changed  by  this  compensation.  The  damping  and  the 
stability  are  both  improved. 

9.  The  addition  of  unity  feedback.  This  makes  a  Type  0  control  out  of  any 
control  type.  Differentiation  added  to  the  feedback  or  the  forward  transfer 
function  of  a  Type  0  control  changes  the  closed-loop  pole  positions.  The 
changed  values  of  the  poles  usually  improve  the  control-system  stability 
and  damping.  If  the  differentiation  is  in  the  forward  transfer  function,  it 
also  adds  a  closed-loop  zero. 

The  incorporation  of  a  feedback  control  system  into  the  towed  body  improves 
the  body’s  dynamics  and  assists  in  the  performance  of  the  towing  missions. 

Body  A  has  been  shown  to  be  dynamically  stable.  It  follows  the  elevator  mid 
cable  tension  step  input  functions  in  longitudinal  motion  with  a  constant  steady- 
state  error.  In  lateral  motion,  the  body  follows  step  input  rudder  or  aileron 
variations  with  constant  error  in  steady  state.  The  towed  body  in  itself  is  a 
Type  0  system  in  both  longitudinal  and  lateral  motion.  For  the  complete  longi¬ 
tudinal  equations  of  motion,  the  transient  solution  consists  of  two  damped 
exponentials  with  very  short  times  to  damp  to  half-amplitude  and  a  damped 
sinusoid  with  a  damping  coefficient  of  0.906.  The  time  to  damp  to  half-amplitude 
of  the  sinusoid  is  0.544  second.  For  the  short,  equations  of  longitudinal  motion, 
the  transient  solution  consists  of  two  of  the  exponentials  of  the  complete 
solution  and  another  damped  exponential  with  a  longer  damping  time.  The  time 
to  damp  to  half-amplitude  of  this  longer-time  term  is  0.866  second.  The  dominant 
or  long-time-to-damp  transient  term  is  largely  determined  by  the  cable  tension 
mid  the  pitch  damping  hydrodynamic  coefficient,  Cm<).  for  both  the  complete 
and  the  short  longitudinal  equations.  The  time  to  damp  can  be  decreased  by 
increasing  me  depression  force  and  decreasing  the  magnitude  of  the  pitch 


1  +  T2s 
1+  Tj's 
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damping  coefficient.  In  the  lateral  equations  of  motion,  the  transient  solution 
has  a  damped  exponential  with  a  very  short  time  to  damp  to  half-amplitude 
(corresponding  to  the  roll  subsidence  mode),  a  damped  sinusoid  with  0.0552 
second  to  damp  to  half-amplitude  (corresponding  to  the  Dutch  Roll  mode),  a 
damped  exponential  with  a  time  to  damp  to  half-amplitude  of  0.2075  second 
(corresponding  to  the  spiral  divergence  mode),  and  a  damped  exponential  with  a 
damping  time  of  1.57  seconds  (cable  mode).  The  term  with  the  largest  damping 
time  is  determined  by  the  cable  tension  and  the  yaw  damping  coefficient,  Cnr 
The  damping  time  can  be  decreased  by  increasing  the  depression  force  and  the 
body  drag,  and  by  decreasing  the  yaw  damping  coefficient.  The  critical  long¬ 
time  damping  terms  for  both  longitudinal  and  lateral  motion  can  be  lowered  by 
increasing  the  depression  force  and  decreasing  the  damping  coefficients, 

Cmq  and  Cnr ■  A  longitudinal  control  system  can  improve  the  damping  and 
stability  and  provide  a  means  of  keeping  the  towed  body  at  a  constant  commanded 
depth.  In  addition,  a  lateral  control  system  can  provide  coordination  so  that 
turning  can  occur  without  skidding  (side  slipping).  The  lateral  control  system 
can  also  be  used  to  orient  the  direction  of  motion. 


LONGITUDINAL  CONTROL  SYSTEMS 

It  has  been  previously  shown  that  Body  A  is  a  Type  0  system  in  that  it 
follows  a  step  input  with  an  error.  The  longitudinal  short  solution  indicates 
that  the  dominant  cable  mode  damps  to  half-amplitude  in  0.866  second.  Here 
we  shall  describe  a  longitudinal  control  system  designed  to  hold  the  towed  body 
at  commanded  depths,  in  addition  to  improving  the  stability  and  damping.  The 
control  system  parameters  depend  on  the  tow  body  with  which  the  system  o|)erntes, 
here  assumed  to  be  Body  A  (which  could  be  improved).  The  effects  of  the 
elements  in  the  control  system  are  given  by  the  transfer  functions  relating  the 
output  to  the  input  of  the  elements  as  functions  of  s,  the  Laplace  transform  of 
time.  The  output/input  conversion  factors  between  angular  and  linear  units  in 
the  elements  is  included  in  the  transfer  function  constants.  The  components  of 
the  control  system  are  assumed  to  be  available.  There  are  many  commercial 
small-sized,  accurate,  and  sensitive  pressure  gauges.  It  is  assumed  that  a  depth 
gauge  converting  depth  to  a  proportional  electrical  signal  is  available.  The 
elevator  servo  can  be  an  electromechanical  or  a  hydraulic  unit,  and  is  represented 
in  the  analysis  by  a  simple  time  lag  and  a  gain  factor.  Idealized  gyroscopes 
without  time  lags  are  assumed.  Commercial  gyroscopes  have  been  developed 
that  are  accurate  and  sensitive  and  are  not  much  affected  by  vibrations,  heat 
variations,  or  pressure  variations.  They  are  available  as  small,  light  units, 
weights  being  about  3  pounds,  lengths  al»ut  5  inches,  and  diameters  about  3 
inches.  The  proposed  control  system  should,  of  course,  be  a  realistic  design. 

Figure  7  consists  of  three  block  diagrams  of  proposed  depth  control 
systems.  The  top  diagram  shows  a  basic  feedback  system  to  control  depth. 

Basic  controls  can  be  adjusted  by  compensating  networks  such  as  minor  loops 
and  filters;  in  the  middle  diagram  a  damping  minor  loop  is  added  to  the  basic 
control,  and  in  the  bottom  diagram  a  filter  is  added  to  the  outer  feedback  path. 
Three  transfer  functions  indicated  in  the  diagrams  are  related  to  the  body  dynam¬ 
ics.  The  transfer  function  relating  the  attitude  angle  rate  with  the  elevator 


angle  deflection,  6/5  e,  is  directly  obtained  from  the  longitudinal  equations  of 
motion.  The  transfer  function  relating  tlie  depth  with  the  attitude  angle  rate, 

A/0,  is  obtained  using  the  vertical  acceleration,  a2: 

JL  is.  J L 

§  =  6e  X  6  \s3 

The  vertical  acceleration  itself  is  calculated  from  the  solution  of  the  iongitudinal 
equations  of  motion: 


where  U  is  the  undisturbed  forward  velocity. 

The  transfer  function  relating  the  depth  with  the  elevator  angle  deflection, 
A/Se,  is  given  by  the  relationship 

_A  _A  _6_ 

6c  =  $  \ 

The  body  dynamics  transfer  functions  for  Body  A,  based  on  the  short-solution 
longitudinal  equations  of  motion,  are  given  by  the  equations: 


Feedback  controls  cnn  be  analyzed  by  many  methods.  The  root  locus 
method  is  convenient  for  repetitive  analysis.  System  design  usually  requires 
several  trial  solutions  that  must  be  optimized.  The  root  locus  plots  indicate 
dumping  and  stability  in  the  location  of  the  roots  of  the  characteristic  equation 
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(closed-loop  poles).  They  show,  in  graphical  form,  variations  of  the  poles  with 
system  gains.  The  zeros  and  the  poles  of  the  open-loop  Iransfer  function  are 
plotted  on  the  complex  plane,  and  root  loci  are  obtained  using  the  transfer 
function  of  the  feedback  control  system: 

C(s)  -  Gls) 

~R(s)  ~  1  +  G(sllUs) 


The  roots  of  the  denominator  are  the  closed-loop  poles.  A  plus  sign  in  the 
denominator  indicates  negative  feedback,  and  a  minus  sign  indicates  positive 
feedback.  Positive  feedback  occurs  when  the  algebraic  sign  of  the  open-loop 
transfer  function  is  a  plus  and  addition  occurs  at  the  error  detector,  or  when  the 
sign  of  G(s)H(s)  is  a  minus  and  subtraction  occurs  at  the  error  detector.  Negative 
feedback  occurs  when  the  product  of  the  sign  of  G(s)IKs)  and  the  sign  of  the 
error  detector  operation  is  a  minus.  The  open-loop  transfer  function  must  be  of 
unit  magnitude  and  of  0-  or  180-degree  phase  at  the  closed-loop  poles.  The 
root  loci  are  plots  of  the  0-degree  phase  requirement  for  positive  feedback,  and 
the  180-degree  requirement  for  negative  feedback.  The  magnitude  requirement 
produces  closed-loop  poles  on  the  root  loci.  A  point  on  the  root  locus  is  a 
closed-loop  pole  when  the  open-loop  transfer  function  is  set  equal  in  magnitude 
to  unity: 


G(s)H(s )  =  1  = 


K(s  -  Z{)  (s  -  Z2) 
(s  -  P{)  (s  -  P2 ) 


The  product  of  the  distances  to  the  open-loop  poles  from  the  point,  divided  by 
the  product  of  the  distances  to  the  open-loop  zeros  from  the  point,  is  equal  to  the 
loop  gain,  R\  corresponding  to  the  closed-loop  pole.  Figure  8A  shows  the  zero 
root  locus  for  the  basic  depth  control  system  corresponding  to  the  upper  block 
diagram  of  figure  7.  Most  of  the  right-hand  locus  is  in  the  right  half  of  the  plane 
and  corresponds  to  instability.  Thu  dominant  closed-loop  pole  for  the  small 
stable  part  of  the  root  locus  has  a  negative  real  part  less  than  0.4.  When  a 
high-pass  filter  is  added  (fig.  8B),  the  root  loci  are  shifted  to  the  left  increasing 
both  the  stability  and  the  damping.  The  dominant  root  now  has  a  negative  real 
part  of  about  4.  The  time  to  damp  to  half-amplitude  decreases  from  about  2 
seconds  for  the  basic  control  alone,  to  about  0.2  second  for  the  basic  control 
plus  the  high-pass  filter.  The  electrical  filter  can  be  designed  using  common 
networks,  for  example: 


O- 

o 


RCs 

1  +  RCs  *  E, 
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E.  MAJOR  LOOP.  WITH  MINOR 
LOOP  AND  FILTER  ADDED 


Figure  S  (Continued). 


Figure  8C  shows  the  root  locus  plot  of  the  minor  loop.  The  dominant  root  is 
determined  by  the  locus  on  the  real  axis  from  -0.8  to  0.  The  value  of  gain 
constant  of  320  produced  a  ciosed-loop  pole  at  —14.21  and  complex  poles  with  a 
damping  constant  of  0.9.  Figure  8D  shows  the  zero  angle  mot  locus  of  the  closet! 
minor  loop  incorporated  into  the  basic  depth  control  system.  The  right-hand 
locus  has  a  very  small  section  in  the  stable  left  half  of  the  complex  plane. 

Figure  8E  shows  the  effect  of  the  addition  of  the  same  high-pass  filter;  the  root 
loci  are  shifted  to  the  left.  The  dominant  pole  can  have  a  negative  real  part  of 
about  7.  The  time  to  damp  to  half-amplitude  can  be  about  0.1  second  with  the 
minor  loop  and  the  high-pass  filter.  Thus,  eonqiensation  can  change  the  control 
system  considerably.  The  transient  response  to  an  input  such  as  a  step  function 
depends  on  the  chosen  compensation.  The  rise  time  to  design  output,  the  over¬ 
shoot  above  the  design  value,  and  the  time  to  damp  are  characteristics  l  aving 
different  damping  requirements.  The  control  system  can  be  optimized  according 
to  several  principles  usually  b-.ed  on  minimizing  errors. 


LATERAL  CONTROL  SYSTEMS 

The  towed  body  is  inherently  stable  due  to  the  stabilizing  effect  of  the 
towing  cable.  Unlike  an  aircraft,  it  tends  to  return  to  its  original  heading  after  a 
disturbance.  The  cable  root  is  dominant  in  the  dynamics.  For  Body  A,  it  has 
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been  show  to  hove  a  half-amplitude  damping  time  of  1.57  seconds.  A  central 
system  can  be  used  to  reduce  the  damping  times.  It  can  be  designed  to  coordinate 
the  rudder  and  the  aileron  controls  so  that  no  side  slip  (skidding)  occurs  diring  a 
turn.  Coordination  can  be  obtained  by  feeding  back  to  the  rudder  input  a  voltage 
proportional  to  the  side-slip  angle  or  the  lateral  acceleration.  Here,  a  rudder 
coordination  computer  is  discussed.  It  may  be  desirable  to  provide  a  method  of 
orienting  the  yaw  angle  of  the  towed  body  by  means  other  than  the  cable.  A 
yaw-rate  body  orientation  control  is  also  discussed. 

The  top  block  diagram  of  figure  9  shows  the  closing  cf  the  rudder  circuit 
through  a  rate  gyroscope  to  form  a  closed  feedback  loop.  The  root  locus  plots 
of  figures  10  (A  through  G)  show  that,  for  the  feedback  ioop,  the  dominant 
closed-loop  pole  is  very  near  the  origin  in  all  cases  for  either  positive  or  negative 
values  of  the  dihedral  coefficient.  When  integration  is  added,  for  example,  by 
use  of  an  integrating  gyro  instead  of  the  rate  gyro,  the  root  locus  plots  indicate 
that  the  closed-loop  poles  can  have  negative  real  parts  considerably  larger  than 
the  dominant  cable  root  particularly  for  the  body  with  the  positive  dihedral 
coefficient. 

The  middle  block  diagram  of  figure  9  shows  connection  of  the  aileron 
circuit  to  the  rudder  circuit  through  a  coupling  network  called  the  rudder  coordina¬ 
tion  computer.  The  method  is  based  on  keeping  the  side-slip  angle  at  zero  during 
turns  by  feeding  back  a  required  rudder  input  for  a  given  aileron  input.'  The  total 
side-slip  angle  is  related  to  the  rudder  and  the  aileron  angles  by  the  equation 

aT  =  rFQ/5Ja  -  TF  (£/8r)r 

where  TF7p/6Q)  is  the  transfer  function  relating  the  side-slip  angle  with  the 
aileron  angle 

TFfp/Sp  is  the  transfer  function  relating  the  side-slip  angle  with  the 
rudder  angle 

When  the  total  side-slip  angle  is  zero,  the  ratio  between  the  perfectly  coordinated 
redder  and  aileron  angles  becomes 

5r_  TFV/5J 

8a  TF  (fi/5r)  rcc 

The  middle  block  diagram  of  figure  9  shows  that  the  negative  of  the  transfer 
function  ratios  is  equal  to  the  transfer  function  of  the  rudder  coordination  com¬ 
puter,  TFrcc.  Unlike  an  aircraft,  positive  redder  requires  positive  aileron,  and 
a  phase  reversal  is  not  required  in  the  rudder  computer  circuit.  The  rudder- 
coordination-computer  transfer  functions  for  Body  A(Cgp  ~)  and  Body  A’  (Cgp  +) 
are  given  by  the  equations 


Body  A 
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YAW  FEEDBACK  CONTROL 


YAW  RATE 
GRYO 


TURN  COORDINATOR 


I  »  6r  I  I  > 

RUDDER  COORDINATION  RUDDER 


COMPUTER  SERVO 


YAW  ORIENTATIONAL  CONTROL  SYSTEM 


Figure  9.  Lateral  control  systems. 
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Body  A' 


Values  of  transfer  functions  that  are  a  decade  in  frequency  away  from  a  reference 
frequency  tire  20  dB  apart  for  first-order  functions.  Thus,  the  rudder-coordination 
transfer  functions  can  be  approximated  for  angular  frequencies  from  about  1  to  10 
radians  per  second  by  the  equations 


Body  A'  (C^ p  +  ) 


The  transfer  function  between  the  total  yaw  rate  and  the  input  voltage  of  the 

rudder-coordination-computer  circuit.  e6  ,  is  given  by  tiie  equation 

a 


=  TF„,  x  TFn-  'Sal  +  TFrs  x  TFrcc  x  TFff  Sr) 

c°a 


where  TFas  is  the  transfer  function  of  the  aileron  servo 
TFrs,  is  the  transfer  function  of  the  rudder  servo 

As  the  equation  above  indicates,  the  two  applied  inputs  of  the  aileron  mid  rudder 
angles  add  algebraically  in  the  yaw-angle  effects.  The  bottom  block  diagram  of 
figure  9  shows  a  yaw  orientation  system  with  rudder  coordination.  An  integrating 
vaw-rate  gyroscope  is  included  The  open-loop  transfer  functions  for  Bodies  A 
and  A'  are  given  by  the  approximate  equations: 


49 


T 
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The  open-loop  transfer  functions  are  based  on  the  approximate  rudder-coordination 
transfer  functions.  The  rudder  produces  the  principal  effects  on  yaw.  The 
orientation  system  can  be  checked  on  analog  computers.  The  more  complete 
rudder-coordination-computer  transfer  functions  can  be  used  if  necessary.  The 
approximate  transfer  functions  can  be  represented  by  the  common  network  below: 


-Wr — - -  - » . . - . O 

R>  k 

[  2  ft2C«+ 1 

-L  IRj  r  R2]Cs  *  * 

_ XL_ 


The  root  locus  plots  for  the  approximate  transfer  functions  are  shown  in 
figure  10  for  the  yaw  orientation  system,  it  can  be  seen,  at  least  at  lower  gains, 
that  both  Bodies  A  and  A’  are  stable.  The  dominant  closed-loop- pole  damping 
time  for  the  positive  dihedral  body  is  about  0.244  second.  For  the  negative 
dihedral  body,  the  damping  time  is  at  the  best  about  0.7  second.  Additional 
networks  can  be  added  to  improve  the  transient  response  if  desired.  Phase-lead 
differentiating  networks  tend  to  narrow  the  system  bandwidth  while  phase-lag 
integrating  networks  tend  to  broaden  the  bandwidth. 

The  control  systems  help  reduce  possible  coupling  between  the  longitudinal 
and  lateral  equations  of  motion.  The  presence  in  the  equations  of  the  two  terms 
listed  below  causes  intercoupiing: 


where  lrr  is  the  body's  moment  of  inertia  about  the  longitudinal  axis 

,  is  the  body's  moment  of  incuiu  about  the  vertical  axis 

/vv  is  the  body’s  moment  of  inertia  about  the  lateral  axis 

The  roll  rate  P  is  not  large  in  the  tow  bodies  due  to  the  cable.  The  addition  of 

control  systems  reduces  the  yaw  rate  R  and  the  pitch  rate  Q. 
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Figure  10. 
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C.  ZERO  ROOT  LOCUS. 

RUDDER  FEEDBACK  LOOP. 
INTEGRATION  ADDED,  C?p~ 
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F.  180°  ROOT  LOCUS. 

RUDDER  FEEDBACK  LOOP, 
INTEGRATION  ADDED.  C.r+ 

Figure  10  (Continued). 
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Figure  10  (Continued). 


In  summary,  in  reference  to  closed-loop  feedback  systems,  it  can  be  said 
that  a  control  system  output  follows  the  input  better,  and  with  less  disturbance 
effects,  and  with  high  forward-transfer-function  gains.  The  output-input  transfer 
function  becomes 

GG. 


1  +  GGf)H 


If  the  gain  G  is  large  and  the  gains  Gh  and  II  are  not  too  small,  the  output-input 
transfer  function  becomes  1  II. 

The  disturbance  output  for  the  network  shown  below  is  given  by 


a, 


1  +  GGfill 
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OUTPUT 


FEEDBACK  NETWORK 


This  becomes  equal  to  the  small  value  1  GH  when  the  gain  G  is  large.  An  input 
filter  called  a  model  can  be  added.  For  a  high  gain,  G,  the  output  follows  the 
model  output.  Adaptive  control  systems  usually  have  high-variable-gain  circuits 
and  still  retain  stability.  They  can  have  electronic  counting  devices  which 
measure  the  gain  and  adjust  to  a  maximum  stable  value.  In  the  case  of  Body  A1 
(Cf(3  +  ).  the  cable  cf  pole  on  the  real  axis  approaches  the  zero  at  very  high 
gains.  The  Dutch  Roll  closed-loop  pole  then  becomes  dominant  at  high  gains  mid 
approaches  the  imaginary  axis;  it  must  be  limited  for  stability. 


TOW  BODY  DESIGN 

Preliminary  design  analysis  has  been  performed  mainly  with  three  bodies, 
designated  Body  A,  Body  I,  and  Body  II.  Detailed  analysis  was  made  in  tioth 
longitudinal  and  lateral  motion  of  Body  A  which  has  a  negative  dihedral  coef¬ 
ficient,  of  Body  A‘  which  has  a  positive  dihedral  coefficient  of  the  same  magni¬ 
tude,  and  of  Body  A"  which  is  of  the  same  form  but  has  a  zero  dihedral  coefficient. 
In  addition  to  the  depressor  angle,  the  wing  location  and  the  tail-ceuter-of- 
pressure  location  in  relation  to  the  body  axis  affect  the  dihedral  coefficient, 
producing  equivalent  dihedral  angles  of  about  .'1  degrees  or  more  depending  on  the 
lateral  area. 

All  the  bodies  studied  have  the  same  .’JtVinch  maximum  diameter,  but  they 
differ  in  other  respects.  Some  important  design  considerations  are  the  drag;  the 
hydrodynamic  coefficients,  particularly  the  damping  coefficients.  Cnr  andCm<;; 
and  the  depression  force.  All  the  bodies  are  axi symmetrical,  streamlined  forms 
with  depressor  and  control  surface  appendages.  We  designed  the  axisymmetrienl 
foims  using  the  David  Taylor  Model  Basin  method  originated  by  Landweber  and 
Gertler.  The  bodies  have  a  length-  diameter  ratio  of  about  2.7  for  minimum  drag. 
Drag  has  an  important  effect  on  the  required  towing  horsepower  and  the  attainable 
depths.  It  has  a  beneficial  effect  on  the  damping  time  of  the  dominant  cable 
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mode  in  lateral  motion.  A  considerable  portion  of  the  required  towing  horsepower 
is  needed  to  counteract  the  body  drag.  For  a  2000-pound  drag,  278  horsepower  is 
required  to  tow  the  body  at  45  knots,  while  for  a  600-foot  towing  cable  with  a 
drag  of  15  pounds  per  foot  1530  horsepower  is  required.  The  attainable  depth  is 
405  feet  at  45  knots  with  a  depression  force  of  14,000  pounds  and  a  2000 pound 
drag. 

The  total  body  drag  depends  considerably  on  the  depressor  and  tail 
surfaces.  Additional  drag  arises  from  the  interference  between  components  and 
from  imperfections  on  the  surfaces.  Hoerner  shows  how  to  estimate  the  total 
drag  on  such  bodies  as  aircraft,  by  adding  together  all  the  drags  of  the  components 
in  drag  areas.  A  drag  area  is  the  product  of  the  drag  coefficient  and  the  area 
that  it  is  defined  for: 


where  q  is  the  dynamic  pressure 
C(jS  is  the  drag  area 

Here,  the  total  drag  is  considered  to  be  the  sum  of  the  depressor  drag,  the  tail 
drag,  the  depressor-induced  drag,  and  the  main  body  drag.  The  imperfection  mid 
interference  drags  are  neglected  since  they  are  small  percentages  of  the  others 
and  can  be  lessened  by  careful  design.  The  drag  area  of  the  induced  drag  is 
approximated  by 


C  2 

_L_ 

ARtx 


where  AR  is  the  aspect  ratio  b~  S 
S  is  the  wing  area 

The  main-body  drag  coefficient  can  be  based  on  wetted  area,  or  cross-sectional 
area  as  it  is  here.  The  drag  coefficient  for  an  axisymmetrica!  body  is  approxi¬ 
mated  by  Hoerner  as  follows: 


where  C ^  is  the  skin-friction  drag  coefficient  of  the*  body 

C’0()  is  the  drag  coefficient  defined  for  the  frontal  area 

The  wing  and  tail  drags  are  based  on  the  plan  area.  The  profile  drag  is  assumed 
to  be  nearly  zero.  The  sectional  drag  coefficient  is  based  on  the  skin  friction 
for  the  upper  and  lower  surfaces  modified  by  factors  that  account  for  the  section 
thickness.  The  coefficient  is  given  by  Hie  equation 
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where  CDs  is  the  sectional  drag  coefficient 
t  is  the  section  thickness 

c  is  the  section  chord 

The  wing  and  tail  thic’-:«ess-to-ch*>rd  ratios  were  taken  to  be  0.2  in  the  calcula¬ 
tions.  The  skin-friction  coefficient  was  assumed  to  he  0.0025.  The  average 
skin-friction  coefficient  is  about  0.<X)25  for  the  body  Reynolds  number  of 
2.5  x  10"  which  corresponds  to  fully  developed  turbulence  over  an  81-inch  body 
at  do  knots  in  50°F  water.  The  transition  to  neutral  turbulence  occurs  on  the 
body  at  about  15  inches  back  of  the  nose.  Transition  to  self-excited  turbulence 
occurs  at  about  22  inches  back  of  the  nose.  Based  on  these  approximations,  tlio 
tow  body  drag  areas  are: 


Main  Body 

13.72  x  Cf*  4.90874 

Drag  Area 

0.1498 

Depressor 

2.48  x  C(  x  6 

0.0372 

Tail 

2.48  x  C{  x  4 

0.0248 

Induced 

0,1188 

0.1188 

Total  Drag  Area 

0.3306 

The  frontal-area  drag  coefficient  is  obtained  by  dividing  the  drag  area  by  the 
frontal  area,  4.90874  square  feet.  The  frontal-area  coefficient  for  the  values 
listed  above  is  0.0673.  The  body  drag  at  45  knots  equals  1910  pounds.  The 
skin-drag  coefficient  may  be  less  than  the  chosen  average  value  of  0.0025  under 
ideal  conditions,  but  the  interference  and  imperfection  drags  and  the  drag  at  the 
tie  point  have  been  omitted.  The  body  is  in  fully  developed  turbulence  at  points 
further  back  than  about  25  inches  from  the  nose  at  45  knots.  For  fully  develop'd 
turbulence,  the  broadband  pressure  at  the  boundary  is  about  80  dB  above 
1  dvne  cm2,  it  is  desirable  to  build  the  body  surface  of  soft  rubber  from  about 
10  inches  back  of  the  nose,  to  dela  turbulence  to  the  self-excited  value.  The 
outer  portion  of  the  body  should  be  made  of  a  low-elastic-moduliis  material  (such 
as  urethane  foam)  of  sufficient  structural  strength,  in  order  to  present  a  filter  to 
the  How-noise  pressure.  Flow  noise  is  present  over  bodies  of  the  tow  lxxl.\  size 
at  45  knots.  It  can  lie  minimized  by  delaying  the  points  of  transition  and  by 
energy  absorption. 

The  depressor  wing  must  be  large  enough  to  produce  a  large  depression 
load  at  low  lift  coefficients.  A  depression  force  of  nearly  I  I. 000  pounds  is 
required  to  attain  depths  of  about  400  feet  with  a  body  drag  of  2000  pounds  and 
a  500-foot  cable  with  a  drag  of  15  pounds  per  foot.  Because  of  induced  drag, 
the  lift  coefficient  should  be  limited  to  values  of  about  0. 1.  The  aspect  ratio 
should  lie  large  for  low  induced  drag.  The  wing  and  tail  section  should  be  thin 
for  small  drag,  but  must  also  hi'  structurally  strong  enough.  Symmetrical  wing 
and  tail  sections  with  a  thickness-to-chord  ratio  limited  to  alxiut  0.2  should  lie 
used.  The  NACA  four  digit  numbering  “vsteni  number  for  this  symmetric  at 
section  is  0020. 

The  horizontal  and  vertical  tail  surfaces  determine  the  damping  coefficients, 
and  C„r.  These  coefficients  affect  all  the  roots  for  both  longitudinal  and 
lateral  rot  ion  fmorabh.  except  the  dominant  entile  roots  for  both  longitudinal 
and  lateral  motion  The  \ allies  arc*  proportional  to  the  square  of  tlu*  distance  from 


the  c.g.  to  the  surface  quarter-chord  point,  and  proportional  to  the  surface  area 
of  either  the  horizontal  or  the  vertical  tail  areas.  The  damping  coefficients 
should  be  large  but  not  so  large  as  to  make  the  cable  roots  too  small. 

The  dihedral  angle  affects  the  roots  of  the  characteristic  equation  and 
hence  the  transient  response.  Positive  dihedral  coefficients  adversely  affect  the 
dominant  cable  root.  The  response  in  yaw  control  systems  appears  to  be  better 
for  positive  dihedral  coefficients.  A  negative  dihedral  coefficient  is  required  to 
pull  up  the  low  wing  in  lateral  gusts  if  the  tie  point  is  connected  at  the  c.g.  of 
the  body.  If  the  tie  point  occurs  above  the  c.g.,  no  dihedral  is  required  due  to 
the-stabiiizing  effect  of  the  towing  cable.  The  diagrams  below  indicate  the  top 
view  and  front  view  of  a  tow  body  with  positive  side  slip  and  a  negative 
dihedral  coefficient. 


The  control  surfaces  are  the  standard  aircraft  type:  an  elevator,  ailerons, 
and  a  rudder.  A  down  elevator  is  assumed  positive,  and  it  produces  a  negative 
moment.  Figure  11  shows  the  effectiveness  factors  of  the  control  surfaces.  A 
0.5  effectiveness  factor, dat/d  5e,  is  provided  by  an  elevator  area  of  about  0.35 
for  an  aspect  ratio  of  7.  The  factor  depends  somewhat  on  the  gap  and  the  type  of 
seal  in  addition  to  the  aspect  ratio,  in  general,  the  variations  are  not  large.  The 
lateral  controls  are  the  ailerons  and  the  rudder.  Positive  aileron,  defined  as 
motion  that  lowers  the  right  wing,  corresponds  to  right  aileron  up.  The  ailerons 
operate  differentially  in  that  when  the  right  aileron  is  up,  the  left  aileron  is  down. 
Positive  rudder,  defined  as  rudder  to  the  left,  produces  a  force  in  the  positive 
Y  direction.  The  ailerons  extend  over  0.7  of  the  depressor  wing  span,  the  ratio 
of  the  chord  of  the  aileron  to  the  chord  of  the  wing  being  0.3.  This  produces 
a  value  of  0.6  for  fj>6o,  the  aileron  moment  coefficient.  The  values  of  this 
coefficient  vary  from  zero  to  about  0.7.  The  ratio  of  ruuder  chord  to  vertical- 
tail-surface  chord  is  chosen  to  be  about  0.35  for  a  0.5  effectiveness  factor.  The 
tail  surfaces  should  have  aspect  ratios  of  about  3. 

For  a  left  turn,  the  tow  body  with  positive  or  left  rudder  requires  positive 
ailerons  to  balance  the  centrifugal  force.  Positive  aileron  lowers  the  right  wing 
and  requires  positive  or  left  rudder.  For  Body  A'  (positive  dihedral  coefficient), 
Qgp  the  rolling  moment  coefficient  due  to  rudder,  opposes  the  proper  bank 
angle:  and  C„ 6a,  the  yawing  moment  coefficient  due  to  aileron  angle,  helps  the 
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Figure  11.  Effectiveness  of  control  surfaces. 
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proper  turn.  For  Body  A  (negative  dihedral  coefficient),  Q5r  helps  the  proper 
bank  angle  whi’v.  C„ga  helps  complete  the  tum.  For  a  30(Vfoot  tum  radius  of 
BodyA,  the  centrifugal  force  is  about  1155  pounds  and  requires  a  bank  angle  of 
about  4.73  degrees. 

4n  addition  to  Bodies  A,  A',  and  A",  a  light  body  designated  Body  I  and 
a  heavy  body  designated  Body  II  were  analyzed  for  a  towing  speed  of  45  knots 
ami  a  straight,  level  path.  The  numerical  parameters  of  Body  I  are: 
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The  hydrodynamic  coefficients  for  longitudinal  motion  of  Body  I  are: 


-4 

-2 

(7 

-0.3 

c. 

-3 

/HOI 

C* . 

-l 

ma 

c 

-8 

zo. 

mq 

-0.125 

-0.5 


The  transfer  functions  for  Body  I  are  shown  in  Appendix  C.  For  Body  II, 
the  numerical  parameters  and  the  hydrodynamic  coefficients  are: 
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The  short-solution  longitudinal  equations  of  motion  indicate  the  times  to  damp 
to  half-amplitude  for  Body  I  with  the  tie  point  at  the  c.g.  to  be  8.295,  0.01752, 
and  0.007255  seconds.  With  the  tie  point  1  foot  above  the  c.g.,  the  damping 
times  of  Body  I  become  0.05584,  0.02436,  and  0.007342  seconds.  For  Body  II 
with  the  tie  point  at  the  c.g.,  the  half-amplitude  damping  times  are  18.15, 
0.02128,  and  0.04667  seconds.  (The  damping  times  for  the  short  solution  of 
Body  A  were  0.866,  0.0368.  and  0.0135  seconds.)  These  results  show  the 
significant  influence  of  the  tie  location  on  the  dominant  cable  root.  The  light 
Eody  I  damps  very  rapidly  for  the  dominant  cable  root  because  of  the  low  value 
of  the  coefficient  Cmq ,  with  the  tie  point  above  the  c.g.  The  weight  and  the 
inertia  are  (optimistically)  too  small.  The  moment  of  inertia  of  Body  II  was 
chosen  to  be  very  large.  This  affects  the  values  of  the  nondominant  roots, 
increasing  their  value  to  about  double  those  of  Body  A.  Figure  12  shows  the 
analog  computer  solution  setup  for  the  longitudinal  equations  of  motion  for 
Body  II. 
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TOW  CABLE 


Tow  cables  have  been  Investigated  by  many,  particularly  M.C.  Eames, 
Figures  13  through  15,  showing  cable  characteristics,  are  from  a  previous  study 
by  the  author.  Figure  13  indicates  that,  at  a  ratio  of  thickness  to  chord  of  0.25, 


t)  2  4  €  8  10  12  14  16  18 

CHORD  THICKNESS  RATIO 


Figure  13.  Cable  cross-section  (frontal)  area  coefficient  versus  chord  thickness 
ratio,  fully  turbulent  flow. 


TION  NUMBER 


the  cable  drag  coefficient  based  on  the  frontal  area  is  0.033.  Table  3 
shows  the  cable  drag  per  foot  at  45  knots  for  three  drag  coefficients  and  three 
cable  thicknesses. 


TABLE  3.  FRONTAL  AREA  DRAG  OF  CABLES 


Cable 

Thickness  (Inches) 

Cable  Drag  for  a 

Frontal  Area  Drag  Coefficient  of: 

0.033  0.035  0.040 

0.375 

5.51 

5.84 

6.68 

0.600 

9.53 

10.11 

11.56 

0.700 

1 

11.12 

12.79 

13.48 

The  drag  coefficients  are  based  on  a  Reynolds  number  of  106  at  normal  incidence. 
The  drag  under  the  varying  conditions  of  actual  motion  will  not  attain  the  values 
shown  in  the  table.  For  a  0.0-inch  thickness  of  optimum  shape,  the  drag  per 
foot  under  operating  conditions  may  be  11  pounds  per  foot  or  more.  A  0.375-inch 
thick,  optimum-shape  cable  with  dynamic  stabilizers  under  operating  conditions 
probably  produces  a  drag  of  8  pounds  per  foot  or  more. 

Figure  16  (A  through  H)  shows  the  body  depth  plotted  against  depression 
force  for  three  depression  force/drag  ratios,  for  cable  lengths  from  500  to  900 
feet,  and  for  cable  drags  of  7  and  15  jxrunds  per  foot.  For  a  cable  drag  of  about 
11  pounds  per  foot,  a  depth  exceeding  400  feet  is  attainable  with  a  (iOO-foot 
cable  length  at  45  knots  for  a  depression  force  of  14,000  pounds  and  a  2000- 
pound  body  drag.  The  attainable  deptn  increases  with  the  depression  force,  but 
the  required  horsepower  to  pull  the  body  increases  with  the  depression  force 
(fig.  17).  About  1560  horsepower  is  required  to  pull  a  body  with  2000-pound 
drag  at  45  knots  if  one  uses  a  (iOO-foot  cable  with  a  drag  of  1 1  pounds  i>er  foot. 
Laminar  flow  cables  are  required  for  minimum  drag  and  cavitation-free  operation 
at  ull  depths  at  this  speed. 
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DEPRESSION  FORCE  IN  THOUSANDS  OF  POUNDS 


Figure  16  (Continued). 
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Figure  16  (Continued). 
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Figure  17.  Required  horsepower  versus  cuble  depression  force;  600-foot  cubic, 
2000-pound  body  drug,  45-knot  speed. 


NOISE  CONSIDERATIONS 

The  tow  body  design  must  consider  the  random,  statistical-type  distur¬ 
bances  that  are  generated  by  surface  activity  and  by  the  turbulent  tow-body 
boundary.  The  recommended  tow  cable  design  feutures  small,  streamlined  cross 
sections.  The  random  energy  generated  in  such  a  cable  should  be  small,  since 
the  boundary  layer  is  laminar  and  no  cavitation  should  occur. 

Most  of  the  tow  body  surface  has  a  turbulent  boundary  at  45  knots.  The 
point  of  neutral  transition  to  turbulence  occurs  at  about  15  inches  back  of  the 
nose,  while  the  self-excited  transition  point  occurs  at  about  22  inches  back  of 
the  nose.  The  value  of  the  broadband  pressure  for  fully  developed  turbulence 
at  the  surface  can  be  estimated.  The  ratio  of  the  broadband  pressure  to  the 
dynamic  pressure  of  the  flow  is  approximately  constant  at  0.0035  in  most  flow 
conditions  in  both  air  and  water.  The  power  spectrum  of  this  flow-noise  energy 
peaks  at  a  particular  value  of  the  Strouhal  number  for  both  air  and  water  at 
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POWER  SPECTRUM, 


normal  speeds.  The  broadband  flow-noise  pressure  expected  at  45  knots  is 
79.72  dB  above  1  dyne  per  square  centimeter.  Figure  17  shows  a  flow-noise 
power  spectrum  curve  obtained  from  pipe  flow.  The  nondimensional  power  spec¬ 
trum  is  plotted  against  the  nondimensional  Strouhal  number.  The  Strouhal  number 
used  is  given  by 


where  V  is  the  flow  velocity 

6  is  the  boundary  layer  thickness 
«  is  the  angular  frequency  in  radians  per  second 

The  power  spectrum  peaks  at  an  angular  frequency  of  5618  radians  per  second 
or  894.12  cycles  per  second  at  45  knots  for  a  0.5-inch  boundary  layer  thickness. 

The  boundary  layer  thickness  can  be  approximated  by  the  nomograms  of 
figures  19  and  20.  Figure  19  gives  the  boundary  layer  thickness  for  both  laminar 
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Figure  18.  Flow  noise  power  spectrum  versus  Strouhal  number. 
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Figure  19.  Nomogram  giving  Reynolds  number  in  water  at  normal  temperature 
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EQUATIONS 
6 . 5.2  x  Rx-V4 
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Figure  20.  Nomogram  giving  smooth,  nut-plate,  boundary  layer  thickness,  for  both  laminar 
and  turbulent  flow 


REYNOLDS  NUMBER  (TURBULENT  FLOW) 


and  turbulent  How  over  a  flat  plate  corresponding  to  the  Reynolds  number  at  the 
point;  the  Reynolds  number  can  be  obtained  from  figure  18.  As  an  example, 
assume  the  distance  is  10  inches  at  45  knots  for  50T  water.  The  Reynolds 
number  becomes  4.4  x  106.  The  boundary  layer  thickness  is  0.17  inch.  The 
peak  of  the  power  spectrum  should  occur  at  16,523  radians  per  second  or  2629.8 
cycles  per  second.  A  representative  boundary  layer  thickness  for  Body  A  is 
about  0.5  to  0.7  inch. 

Random  energy  generated  at  the  surface  by  ship’s  motion  and  wave  motion 
applies  an  input  to  the  towing  cable.  The  spectrum  of  this  energy  is  probably 
fairly  wide-band.  The  toy  ing  cable  conducts  some  of  this  energy  down  the 
cable  in  the  form  of  compressional  and  transverse  waves,  the  cable  acting  as  a 
filter.  The  average  input  power  to  the  cable  is  given  by 

1  f  «o 

P ;  -  —  /  Q>(oi)da> 

1  2tc 

where  <D(co)  is  the  input  power  spectrum. 

'1  he  power  output  from  the  cable  is  g>’  -n  by 

Pq  -  /  K»c|  2  Ofcol  du> 

where  Gc  is  the  transfer  function  of  the  cable. 

If  the  cable  is  represented  by  a  simple  low-pass  filter  with  a  transfer  function  of 

1 

1  +  RCs 

and  the  input  power  spectrum  is  white  noise,  the  output  power  becomes 

K 

2  RC 

where  K  is  the  constant  white-noise  spectrum. 

The  towing  cubles  probably  filter  out  much  of  the  input  surface  energy 
even  at  frequencies  below  1000  cycles  per  second. 


RESULTS  AND  CONCLUSIONS 

The  following  results  and  conclusions  are  based  on  a  mathematical 
analysis  of  the  motion  of  underwater  towed  bodies: 

1.  A  practical  tow  body,  designated  Body  A,  was  designed.  The  maximum 
diameter  is  ,30  inches ,  and  the  main-body  length/diameter  ratio  is  2.7. 

2.  The  tie-point  location  was  shown  to  be  important  to  the  body’s  stability. 
It  should  bo  forward  and  above  the  body’s  center  of  gravity. 
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3.  The  hydrodynamic  coefficients  can  be  predetermined  to  a  fair  degree  of 
accuracy.  A  mathematical  model  can  be  set  up  which  can  be  checked  on  analog 
or  digital  computers. 

4.  The  most  important  parameters  determining  the  dominant,  characteristic 
equation  roots  for  both  longitudinal  and  lateral  motion  are  the  cable  tension  and 
the  damping  coefficients,  Cmq  and  Cnr.  The  root  values  determine  the  body’s 
stability  and  transient  response. 

5.  The  body  drag  is  an  important  parameter  determining  the  required  towing 
horsepower  and  the  attainable  depth.  The  total  body  drag  is  the  sum  of  the  drags 
due  to  the  main  body,  the  depressor  wing,  the  tail  surfaces,  and  induced  drag. 
Induced  drag  is  a  considerable  portion  of  the  total  drag.  The  minimum  attainable 
drag  under  idealized  conditions  is  probably  about  1750  pounds  at  45  knots, 

6.  Tow  bodies  can  be  easily  designed  to  have  longitudinal  and  lateral 
dynamic  stability.  The  hydrodynamic  coefficients  have  the  major  influence  on 
the  roots  of  the  characteristic  equations. 

7.  The  tow  bodies  without  automatic  controls  are  “Type  0”  systems.  They 
follow  step  inputs  of  the  control  surfaces  with  a  following  error. 

8.  The  depressor  wing  area  should  be  sufficient  to  produce  high  depression 
forces  at  low  value  of  lift  coefficients.  Low  lift  coefficients  and  high  aspect 
ratios  reduce  the  induced  drag. 

9.  A  depressor-wing  dihedral  angle  is  required  if  the  tie  point  is  at  the  body’s 
center  of  gravity.  A  negative  dihedral  coefficient  is  more  favorable  than  a 
positive  coefficient  in  its  effect  on  the  dominant  cable  root.  A  better  lateral 
control  system  can  be  obtained  with  a  positive  coefficient. 

10.  In  choosing  the  tow  body  construction  materials,  one  should  consider  the 
pressure  fluctuations  of  random  flow  noise. 

11.  The  control  surfaces  can  be  designed  with  sufficient  area  to  provide  an 
acceptable  degree  of  control  effectiveness. 

12.  Control  systems  greatly  improve  the  transient  characteristics.  The  stability 
and  damping  can  be  favorably  changed. 

13.  A  longitudinal  control  system  is  presented.  This  system,  in  addition  to 
holding  the  tow  body  at  a  commanded  depth,  improves  the  damping  and  stability. 

14.  A  lateral  control  system  is  presented.  This  system  improves  the  stability 
and  damping,  coordinates  the  rudder  and  ailerons  so  that,  no  skidding  occurs 
during  turns,  and  provides  a  yaw-rate  directional  control  for  the  tow  body. 

15.  Intercoupling  of  the  longitudinal  mid  lateral  motion  is  very  small  even 
without  automatic  controls.  Roll  damping  of  the  tow  body  itself  is  excellent. 

16.  Turning  is  considerably  aided  by  the  introduction  of  turn  coordination. 

The  required  bank  angle  is  less  than  5  degrees  for  a  300- loot-rad > us  turn  for 
Body  A  at  45  knots. 

17.  Flow  noise  study  indicates  a  broadband  pressure  of  80  dB  above  1  dyne 
per  cm2  in  the  turbulent  boundary  layer  of  the  tow  body.  A  nondimensional  power 
spectrum  curve  is  presented.  The  frequency  at  which  the  power  spectrum  peaks 
depends  on  the  speed  and  the  boundary  layer  thickness.  Nomograms  are 
presented  which  can  be  used  to  estimate  the  boundary  layer  thickness. 

18.  The  streamlined  towing  cable  attenuates  much  of  the  higher-frequency 
random  energy  travelling  down  the  cable,  particularly  for  transverse  modes. 

19.  The  attainable  depths  and  required  horsepower  depend  primarily  os. 
towing  cable.  The  tow  body’s  depression  force/drag  ratio  must  be  large. 
Small-cross-section,  laminar-flow  streamlined  cables  are  required  for  the  low 
drag  per  foot,  and  for  the  prevention  of  cavitation  at  all  depths. 
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Figure  A-l.  Angle  of  attack  -  elevator  angle;  short  solution. 
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figure  A  2  Attitude  angle  _  elevator  angle;  short  solution. 
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Figure  A-5.  Attitude  angle  -  elevator  angle;  complete  longitudinal  solution. 
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Figure  A-ti  Angle  of  attack  -  elevator  angle;  complete  longitudinal  solution. 


Figure  A-H  Angle  of  nltnck  -  cubic  tension;  complete  longitudinal  solution 
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Figure  A- 11.  Side-slip  angle  -  rudder  angle. 
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Figure  A- 12.  Yaw  angle  -  rudder  angle. 
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